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1.0  Simmiary. 


Biodynamic  Research  Corporation  (BRC)  completed  an  SBIR  Phase  I  project  to  study  the 
feasibility  of  developing  a  model  of  the  Aviator’s  Breathing  System  (ABS).  The  motivation 
for  the  project  was  the  desire  to  develop  a  model  which  could  simulate  the  cardiovascular 
and  respiratory  responses  to  altitude  and  acceleration  stress  encountered  in  high  performance 
military  aircraft.  Once  validated,  the  model  could  be  employed  to  integrate  existing 
physiologic  knowledge  and  to  estimate  the  physiologic  effects  of  changes  in  the  environment 
or  the  design  of  aircrew  protective  equipment. 

The  Phase  I  effort  was  concentrated  in  developing  computer  software  modules  which  were 
capable  of  simulating  the  main  effects  of  environmental  variables  on  the  breathing  ^stem 
hardware  and  respiratory  and  cardiovascular  systems.  Software  modules  were  developed  and 
tested  for  simulation  of:  (1)  the  flows  and  pressures  within  the  breathing  gas  delivery 
system;  (2)  the  flows,  pressures,  and  respiratory  gas  distribution  within  the  lung;  and  (3)  the 
steady-state  flows  and  pressures  within  the  cardiovascular  system.  Subprograms  were  also 
developed  to  compute  altitude-barometric  pressure  relationships  as  well  as  passenger  cabin 
pressures  in  military  aircraft. 

In  addition  to  the  software  development,  BRC  reviewed  and  organized  the  Government 
furnished  data  from  a  series  of  manned  rapid  decompressions  known  as  the  “EONS 
Experiments.  ”  The  data  from  approximately  170  experimental  decompressions  were 
screened  for  their  suitability  for  use  in  parameter  selection  and  validation  of  the  respiratory 
modelling  software.  The  data  appears  to  be  highly  coherent  and  fully  usable  for  model 
validation.  The  summarized  results  of  the  EONS  Experiments  represent  new  physiologic 
database  which  is  directly  applicable  to  USAF  operational  scenarios.  The  respiratory  gas 
composition  time  histories  represent  original  scientific  data  suitable  for  publication  and  use 
by  the  aerospace  medicine  community  independently  of  their  value  in  the  ABS  modelling 
effort. 

The  work  accomplished  in  Phase  I  has  laid  the  groundwork  for  developing  an  integrated  ABS 
Model  in  a  Phase  II  SBIR  program.  BRC  has  written  prototype  softw^e  modules  which 
form  the  major  building  blocks  for  the  Phase  II  ABS  Model.  We  conclude  that  the 
development  of  an  integrated  ABS  Model  is  feasible  and  desirable.  We  also  conclude  that 
there  is  suflBcient  existing  data  for  selection  of  parameters  and  independent  validation  of  the 
major  components  of  the  Phase  II  ABS  model.  However,  the  complete  validation  of  the 
combined  effects  of  acceleration,  pressure  breathing  and  altitude  on  the  ABS  must  await  new 
experimental  data.  Fortunately,  the  attempt  to  simulate  the  physiologic  response  to  combined 
stressors  will  explicitly  point  to  missing  data  and  guide  the  development  of  experimental 
protocols  to  elucidate  their  significance. 
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SECTION  2.0 


2.0  Introduction. 


This  final  report  concludes  an  effort  by  Biodynamic  Research  Corporation  (BRC)  of 
San  Antonio,  Texas  to  conduct  a  study  to  investigate  the  feasibility  of  creating  a  computer 
model  of  the  respiratory  gas  dynamics  in  the  military  aviator’s  breathing  system.  The  work 
was  conducted  under  Contract  #F41624-93-C-6009  through  the  Armstrong  Laboratory  and  is 
entitled  “Modelling  Respiratory  Gas  Dynamics  in  the  Aviator’s  Breathing  System." 

The  aim  of  the  study  was  to  establish  the  feasibility  of  creating  a  comprehensive 
mathematical  model  of  the  Aviator’s  Breathing  System  (ABS)  which  could  be  employed  to 
simulate  the  dynamics  of  respiratory  gas  movement  and  exchange  in  breathing  systems  for 
high  performance  military  aircraft.  The  ABS  Model  could  be  employed  as  a  tool  for 
integrating  existing  physiologic  knowledge  and  for  guiding  future  research  by  ocplicitly 
identifying  the  new  data  required  to  more  fully  describe  the  physiologic  consequences  of 
breathing  system  design,  as  well  as  assist  in  elucidating  the  major  effects  of  the  environment 
and  protective  systems  on  aircrew  performance,  safety,  and  comfort.  The  major  technical 
objectives  of  the  Phase  I  project  are  listed  below: 

(1)  Conduct  a  review  of  the  literature  to  assess  previous  models  for  their  suitability  as 
components  of  the  ABS  model  and  to  glean  physiologic  data  from  the  literature  for 
use  in  parameter  identification  and  model  validation. 

(2)  Synthesize  mathematical  models  of  the  components  of  the  Aviator’s  Breathing 
System  including:  pulmonary  mechanics  and  gas  exchange,  cardiopulmonary  gas 
transport,  and  parametric  models  of  oxygen  delivery  system  components. 

(3)  Analysis  of  Government  furnished  experimental  rapid  decompression  data. 

(4)  Develop  model  parameters  and  partially  validate  tentative  models  of  the  ABS 
components, 

(5)  Assess  the  feasibility  of  development  of  a  comprehensive  ABS  model  in  Phase  n. 

Biodynamic  Research  Corporation  performed  all  of  the  analysis  for  this  contract.  This  final 
report  documents  the  work  generally  according  to  the  objectives  listed  above.  Section  3 
provides  general  background  on  modelling  and  simulation  of  relevant  physiologic  systems 
and  breathing  system  components.  Sections  4  through  6  discuss  the  component  system 
models  separately  and  give  additional  background  and  technical  details.  Section  7  presents 
the  analysis  of  the  Government  furnished  experimental  data  which  is  summarized  in 
Appendix  A.  Finally,  Section  8  presents  conclusions  and  recommendations  together  with  the 
assessment  of  the  feasibility  of  creating  a  comprehensive  integrated  ABS  Model  in  Phase  H. 
In  addition  to  the  experimental  data  summary,  the  appendices  to  this  report  contain  a 
comprehensive  bibliography  of  relevant  literature  and  the  listings  of  the  computer  programs 
developed  during  the  project. 
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3.0  Technical  Background  and  Literature  Review. 

3.1  Modelling  Oxygen  Delivery  Systems. 

There  are  ample  laboratory  data  describing  the  flow/pressuie  relationships  within  breathing 
systems. However,  other  than  the  present  project,  the  authors  are  not  aware  of 
any  attempts  to  model  USAF  breathing  gas  delivery  systems.  There  is  one  rqx>rt  that 
describes  a  model  of  a  simple  US  Navy  demand  regulator,  but  no  details  of  the  system 
equations  were  given.*  The  breathing  system  model  created  during  the  present  study 
simulates  the  pressures  and  flows  in  a  generic  oronasal  mask  during  breiuhing  through  a 
diluter-demand  breathir.g  regulator. 

Because  breathing  regulators  are  inherently  nonlinear  devices,  writing  and  solving  models  of 
their  function  is  difficult.  One  alternative  method  of  modelling  regulators  is  to  characterize 
their  performance  parametrically  over  the  operating  range  of  interest.  >^th  the  exception  of 
the  Navy  study  cit^  above,  parametric  modelling  of  empirical  data  has  invariably  b^  used 
for  describing  regulator  performance.  A  more  rigorous,  but  much  more  difficult  method  is 
to  write  and  solve  the  system  equations  for  the  regulator  of  interest.  For  the  Phase  I ABS 
study,  the  regulator  was  modelled  parametrically,  while  the  oronasal  mask  was  modelled  by 
solving  the  system  equations  describing  its  performance. 

3.2  Modelling  Respiratory  Function. 

Models  of  steady  state  respiratory  gas  exchange  are  well  established.’  The  most  commonly 
employed  of  these  is  the  Alveolar  Gas  Equation  (AGE)*'*  which  describes  the  steady  state 
relationship  between  barometric  pressure,  fraction  of  inspired  oxygra,  and  the  fra^ons  of 
oxygen  and  carbon  dioxide  in  the  alveoli.  Although  useful  in  many  aerospace  applications, 
the  AGE  cannot  describe  the  breath-to-breath  transients  which  occur  during  ordinary 
breathing,  nor  can  it  be  employed  to  describe  the  dynamics  of  respiratory  gas  exchange 
during  rapid  changes  in  ambient  pressure  or  when  ffiere  is  a  net  exchange  of  nitrogen 
between  the  lung  and  the  environment.  During  most  flight  operations,  even  at  high  altitudes, 
as  long  as  the  aircraft  cabin  pressurization  system  remains  functional,  the  AGE  may  be 
employed  to  set  the  breathing  system  design  criteria  for  avoiding  undue  hypoxia.  However, 
during  rapid  changes  in  cabin  ^titude  and/or  changes  in  the  pressure  or  composition  of  the 
inspired  breathing  gas,  the  assumptions  required  for  use  of  the  AGE  are  violated.  Thus, 
there  will  be  significant  errors  between  the  actual  alveolar  gas  composition  and  that  predicted 
by  the  AGE  until  the  respiratory  steady  state  is  re-established. 

To  adequately  describe  breath-by-breath  gas  dynamics  in  the  lung  during  pressure/ 
composition  transients,  a  model  which  includes  the  spatial  and  temporal  variation  of 
ventilation  and  perfusion  within  the  lung  is  required.  There  are  several  models  described  in 
the  literature  which  describe  cardiopulmonary  function  in  clinically  significant  pulmonary 
disease  or  for  teaching  respiratory  physiology. Unfortunately,  these 
were  developed  for  modelling  a  particular  physiologic  response  or  for  a  particular 
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Figure  3*1 .  Aviator's  Breathing  System 


environment  or  clinical  situation.  For  the  most  part,  previous  models  are  not  easily 
adaptable  to  the  high  performance  aviation  scenario  because  they  do  not  describe  many  of  the 
relationships  between  the  environment,  the  breathing  system,  and  the  crew  member’s  other 
protective  systems.  As  an  alternative,  BRC  proposed  Ae  development  of  an  ABS  model 
tailored  to  model  the  gas  flows  and  composition  within  the  breathing  and  respiratory  systems 
in  typical  operational  and  experimental  environments  encountered  by  military  aircrews. 
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For  the  ABS  Phase  I  project,  two  modules  were  developed  which  together  describe  the 
vertical  ventilation  pattern  in  the  lungs  and  the  concentration  of  a  tracer  gas  throughout  the 
lungs  during  inspiration  and  expiration.  The  lung  is  modelled  as  four  regions  of  equal  mass 
that  are  located  along  the  verti^  axis  of  the  lung.  The  model  accounts  for  the  volume  and 
ventilation  differences  created  by  gravity,  which  have  traditionally  been  called  ‘‘interregional 
differences”  by  physiologists.  A  detail^  description  of  the  ABS  Phase  I  models  is  given  in 
Section  5.0. 

3.3  Modelling  the  Cardiovascular  System. 

Ultimately,  the  ABS  model  seeks  to  describe  the  cardiopulmonary  physiology  associated  with 
the  various  environmental  stressors  associated  with  high  performance  flying.  Early  in  the 
Phase  I  program,  it  became  apparent  that  a  reasonably  complete  model  of  the  cardiovascular 
system  would  be  required  to  adequately  model  the  effects  of  positive  pressure  breathing 
(PPB)  for  altitude  and  acceleration  protection  (PBA  and  PBG).  Thus,  it  was  decided  to 
create  a  steady  state  cardiovascular  model  for  the  Phase  I  program  which  could  be  extended 
to  a  dynamic  model  in  Phase  II.  We  concentrated  our  literature  search  on  modem 
cardiovascular  modelling  work  which  appeared  to  be  suflSciently  comprehensive,  but  not 
overly  complicated.  In  particular  the  cardiovascular  model  described  in  Section  6.0  parallels 
the  developments  presented  in  a  paper  by  White,  et  al.**  It  also  draws  from  papers  by 
Jaron,  et  al.,'’’^-^‘  whose  approach  was  similar  to  that  presented  in  earlier  papers  by 
Rideout,^  Snyder,^  and  Avula.^  All  of  these  authors  cited  the  earlier  work  of 
Womersley^  who  developed  equations  describing  the  convection  of  fluids  in  elastic  tubes  as 
well  as  those  describing  the  propagation  of  pressure  waves  and  pressure  dependent  changes 
in  tube  radius. 

3.4  Integrated  Models  of  Cardiopulmonary  Function. 

The  authors  are  aware  of  only  one  integrated  model  of  the  physiologic  function  of  the  human 
body,  a  computer  model  known  as  Human^*.  Human  is  a  teaching  program  for  medical 
students  which  is  an  integrated  computer  model  of  the  entire  body.  BRC  has  experience  with 
an  old  version  and  has  been  impressed  with  its  fidelity,  given  its  simplicity.  BRC  has 
purchased  the  source  code  for  Human  and  has  working  versions  written  in  both  MS-DOS 
Fortran  and  MS-DOS  Quick  Basic  (QBX).  The  general  approach  employed  in  Human  is  to 
model  physiologic  function  by  use  of  curve  fits  to  empirical  data.  This  approach  limits 
Human’s  use  to  conditions  which  have  been  incorporated  in  the  model.  Unfortunately, 
Human  as  presently  configured  only  allows  small  excursions  from  sea  level,  + 1  Gz 
conditions.  However,  the  approach  employed  in  Human  to  implement  physiologic  control 
systems,  such  as  the  baroreceptor  system  for  control  of  cardiac  output,  is  simple  and 
effective.  It  is  suggested  that  a  similar  approach  can  be  employed  to  model  the  physiologic 
feedback  control  in  the  Phase  II  ABS  model.  It  would  be  relatively  straightforward  to  refine 
and  modify  the  Human  code  to  simulate  the  range  of  conditions  required  for  the  ABS  model. 
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3.5  Modelling  the  Environment. 

By  the  environment,  we  mean  the  physical  environment  including  such  things  as  the 
barometric  pressure,  the  temperature,  the  acceleration,  etc.  Most  environmental  variables 
would  be  implemented  as  inputs  to  the  hardware  and  physiologic  models.  The  Phase  I  ABS 
Model  includes  subroutines  which  model  the  variation  in  atmospheric  pressure  with  altitude 
as  well  as  the  change  in  MILSPEC  cabin  pressure  with  altitude.  Both  models  are  included  as 
part  of  the  model  of  the  regulator  and  oronasal  mask  described  in  Section  4.0. 

3.6  Summary. 

At  the  present  time,  there  are  no  comprehensive  physiologic  models  which  can  describe  the 
physiologic  consequences  of  changes  in  operating  conditions  and  design  of  the  Life  Support 
Systems  (LSS)  employed  in  modern  high  performance  aircraft.  BRC  has  created  tentative 
computer  models  of  the  components  required  to  implement  an  integrated  model  of  the  ABS 
which  simulates  the  function  of  the  hardware  and  the  physiologic  subsystems  in  response  to 
environmental  changes.  Detailed  description  of  the  newly  developed  models  are  given  below 
with  comments  on  the  improvements  necessary  to  create  an  integrated  ABS  Model  in  a  Phase 
II  SBIR  program. 
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4.0  Mask  and  Regulator  Models. 

4.1  General  Description. 

The  mask  and  regulator  models  simulate  the  time  varying  flows  and  pressures  within 
breathing  ^stem  from  the  regulator  outlet  to  the  trachea  during  the  breathing  cycle.  Both 
models  are  based  on  the  performance  standards  promulgated  by  the  Air  Standardization  and 
Coordinating  Committee  in  Air  Standard  61/21.” 

4.2  Mask  Model. 

The  mask  is  modelled  as  a  dead  space  and  two  variable  area  orifices  through  which 
inspiratory  and  expiratory  fiows  separately  pass.  The  inspiratory  valve  model  relates  the 
flow  and  pressures  between  the  mask  supply  hose  and  the  oronasal  cavity  during  inspiration. 
The  expiratory  valve  model  relates  the  flow  and  pressure  drop  between  Ae  mask  cavity  and 
the  external  ambient  atmosphere  during  expiratory  flow.  Both  valve  models  are  based  on 
empirical  data  collected  on  the  RAF  P/Q  mask,^^  which  meets  the  Air  Standardization 
Coordinating  Committee  flow  resistance  standards.”  Other  valve  models  can  be  substituted. 
The  general  form  of  the  model  is 

f'  -  ■  [“-n 


where,  is  the  mask  valve  area,  is  the  instantaneous  volume  flow  through  the 
valve.  The  model  actually  calculates  the  mass  flow  through  the  valves  and  converts  the 
result  to  volume  flow  using  the  gas  composition  and  the  prevailing  temperature  and  absolute 
pressure.  The  valve  area,  A^,  is  in  turn  a  function  of  the  pressure  drop,  across  the 

valve.  The  valve  models  are  solved  by  estimating  the  area  of  the  valve  from  the  current 
pressure  conditions  and  then  computing  the  flow  through  the  valve  using  ideal  orifice 
equations.  Mask  leaks  can  be  modelled  by  parallel  fiow  paths  to  ambient,  but  none  are 
included  in  the  present  model.  Expiratory  valve  compensation  is  simulated  by  adding  any 
positive  difference  between  mask  hose  pressure  and  ambient  pressure  to  the  expiratory  valve 
down  stream  pressure.  In  the  absence  of  compensation,  the  down  stream  pressure  for  the 
expiratory  valve  is  ambient  atmospheric  absolute  pressure  plus  the  valve  cracking  pressure. 

The  mask  hose  and  connectors  are  modelled  as  simple  tubes  with  flow-pressure  drop 
relationships  based  on  quadratic  fits  to  typical  data.  Data  describing  specific  actual 
equipment  can  be  substituted. 

4.3  Regulator  Model. 

The  authors  know  of  no  published  functional  models  of  modem  breathing  regulators.  There 
is  one  report  that  describes  a  model  of  a  simple  US  Navy  demand  regulator,  but  no  details  of 
the  system  equations  were  given.*  Transfer  function  models  of  breathing  regulators  have 
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also  been  employed  to  investigate  regulator  stability.’  However,  breathing  regulators  are 
inherently  nonlinear  devices,  so  that  a  transfer  function  approach  can  only  accurately  model 
regulator  function  over  a  limited  range  of  operating  conditions.  Therefore,  we  chose  to 
simulate  the  performance  of  an  ASCC  regulator  rather  than  its  function.  The  regulator  is 
simulated  parametrically  as  a  mass  flow  into  the  distal  end  of  the  mask  tube  in  response  to 
pressure  changes  at  the  regulator  outlet.  Pressure  Breathing  for  Altitude  (PBA)  or  for  high 
acceleration  (PBG)  are  accommodated  by  a  regulator  outlet  reference  pressure  which  can  be 
made  a  function  of  ambient  pressure  or  acceleration.  The  outlet  gas  composition  may  be 
included  in  the  simulation  to  duplicate  the  regulator  dilution  performance  for  various 
environmental  conditions.  The  Phase  II  regulator  model(s)  can  be  modified  to  rq>resent  the 
specific  performance  of  any  breathing  regulator  for  which  empirical  flow-pressure  data  are 
available. 

4.4  Airway  Resistances. 

In  the  mask  model,  airway  resistances  are  modeled  by  simple  quadratic  fits  to  physiologic 
data.  The  resistance  of  the  bronchi  and  trachea  are  included  in  the  model,  as  is  the  pressure 
loss  in  the  oronasal  cavity.  The  intrapulmonary  pressures  estimated  in  the  Mask  Model 
represent  the  mean  pressure  at  the  entrance  to  the  alveolar  space.  The  Pulmonary  models 
described  in  Section  5.0  compute  the  regional  intrapulmonary  pressures  more  accurately  than 
the  technique  employed  here.  The  Phase  II  ABS  model  will  integrate  the  mask  and 
pulmonary  models  and  will  employ  the  more  accurate  methods  described  in  Section  5.0,  but 
simple  models  were  included  in  the  Phase  I  Mask  Model  in  the  interest  of  completeness. 

In  the  Mask  Model,  the  resistances  for  mouth  breathing  and  nose  breathing  are  accounted  for 
separately.  The  following  empirical  model”  is  employed. 


AP  =  ATj-K  +  .  [4-2] 

The  values  for  the  K’s  were  also  taken  from  reference  29.  For  mouth  breathing, 

Kj  =  2.35-10*  Pa-sec-m*  and  Kj  =  2.94-10’  Pa-sec-m"*. 

For  the  nose  breathing  there  are  separate  K2’s  for  inspiration  and  expiration,  but  K,  is  the 
same  for  both  phases. 


K,  =  2.94-10*  Pa-sec-m*,  and 
Inspiratory  Kj  =  2.94-10*  Pa-sec-m'*,  and 
Expiratory  Kj  =  3.92-10*  Pa-sec-m"*. 
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The  overall  pressure  drop  is  modelled  by  assigning  the  relative  fraction  of  the  total  breathing 
flow  to  the  mouth  and  nose,  respectively.  The  individual  drops  are  weighted  by  the  flow 
fractions.  Thus, 

F.  ^  F,  •  1  ,  14-3] 

and 


•AP 


+  F.AP. 


[4-4] 


4.5  Forcing  Function. 

The  Phase  n  ABS  Model  will  employ  the  Pulmonary  Model  to  goierate  mask  flows  in 
response  to  diflferences  in  intrapleural  pressure  and  external  pressures  on  the  thorax  and 
abdomen.  For  the  present,  we  employed  a  simple  sinusoidal  forcing  function  for  the  Mask 
Model  to  simulate  the  volume  change  generated  by  the  lung  during  breathing.  The  presrat 
coding  models  the  lung  volume  change  by  the  following  fimction. 


y 

»i(»)  =  MF^  -  -^•Cos(<.t)  • 

SO  that  the  lung  flow  is 

Vj^  =  .  [4-6] 


4.6  Simulation  Results. 

To  demonstrate  the  model  output,  system  flows  and  pressures  were  computed  for  moderate 
breathing  at  ground  level.  The  respiration  rate  was  20  breaths'min  *,  the  tidal  volume  was 
1.5  L,  and  the  peak  inspiratory  flow  was  94.2  L-min  ’.  Figure  4-1  shows  the  system  flows 
for  the  regulator,  mask,  and  lung.  Positive  flows  are  into  the  mask  and  lung  (and  out  of  the 
regulator),  while  n^ative  flows  are  flows  out  of  the  mask  and  lung.  An  inter-breath  pause  is 
simulated  by  a  short  period  between  cycles  during  which  there  is  no  flow.  The  flows  in  the 
mask  and  lung  are  computed  at  Body  Temperature  and  Pressure,  Saturated  (BTPS) 
conditions,  while  the  flow  from  the  regulator  outle  trough  the  mask  tube  are  computed  at 
Ambient  Temperature  and  Pressure,  Dry  (ATPD)  conditions.  The  mask  and  lung  flows  are 
virtually  identical,  so  that  their  curves  are  superimposed  in  the  plot.  Figure  4-2  illustrates 
the  system  pressures  on  both  absolute  and  differential  pressure  scales.  The  pressures 
correspond  to  the  flows  shown  in  Figure  4-1.  Thus,  the  minimum  pressures  correspond  to 
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.  System  Rows  for  the  Reguletor,  Mask,  and  Lung 


[U|UI/«J«||-|]  MO|d 
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.  System  Pressures  on  Absolute  end  Differentiel  Pressure  Sceles 
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the  peak  inspiratory  flows,  while  the  maximum  pressures  corre^xmd  to  the  peak  expiratory 
flows.  Figure  4-3  shows  the  system  volumes  (cumulative  flows).  As  shown,  the  flows  into 
and  out  of  the  mask  and  lung  are  equal  for  a  single  cycle  of  the  sinusoidal  breathing  pattern. 
Thus,  the  mask  and  lung  volumes  rise  and  fall  around  their  initial  values.  Since  the  r^ulator 
flow  is  unidirectional,  the  volume  increases  during  each  inspiratory  cycle  and  remains 
constant  during  the  expiratory  cycles. 

The  model’s  ability  to  simulate  a  rapid  decompression  is  illustrated  in  Figures  4-4,  4-S,  and 
4-6.  Figure  4-4  displays  the  system  absolute  pressures  during  a  S  psig  (3.4S  kPa) 
decompression  from  sea  level  to  11,(X)0  feet  (3353  m)  altitude.  The  diffemitial  pressures  are 
given  in  Figure  4-5.  The  “breathing”  in  the  model  was  paused  during  the  1.65  sec  rapid 
decompression.  The  model  contains  a  simulated  relief  v^ve  which  controls  the  r^ulator 
outlet  and  mask  tube  pressure.  As  implemented,  the  relief  is  ideal.  That  is,  it  has  no  flow 
resistance.  Thus,  the  mask  cavity  pressure  rise  is  limited  during  the  simulated 
decompression.  The  Phase  n  regulator  model  will  include  code  to  simulate  the  function  of 
an  actual  compensated  poppit  relief  which  will  give  a  more  realistic  response  to  pressure 
changes.  Figure  4-6  shows  the  system  flows  corresponding  to  the  pressure  transients  in 
Figure  4-5. 
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Figur*  4-3.  Cumulativ*  Voluin* 


ID 


•uiniOA  •Aininiuno 

^  M  V-  e 


%  I  % 

e  e  e 

[CyUi]  •uinfOA  •ARtinuino 


110000 


[Oh  uiui]  tjntMJd  tiniotqv 


40 


41 


[3M]  •UIIJ. 


System  Flows 


[U|UI/«J«}ii]  AU>y 


Time  [sec] 


SECTION  5.0 
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5.0  Pkilmooary  Model. 

5.1  Overview. 

The  Pulmonary  Model  describes  the  vertical  ventilation  pattern  in  the  lungs  and  the 
concentration  of  a  tracer  gas  throughout  the  lungs  during  a  breath.  Two  computerized 
models  comprise  the  Pulmonary  Model.  In  these  models,  each  lung  is  made  up  of  four 
regions  of  equal  mass  that  are  located  at  different  vertical  positions.  Figure  S-1  gives  an 
overview  of  the  inputs  and  ouq)uts  of  the  two  models.  The  first,  the  Vratilation  Model, 
determines  the  ventilation  pattern  along  the  vertical  axis  of  the  lung.  The  results  of  this 
model  are  input  into  a  second  model,  the  Convection-Diffusion  Model,  which  calculates  the 
concentration  of  a  tracer  gas  throughout  each  lung  during  a  breath.  The  Pulmcmary  Model  is 
designed  to  simulate  experiments.  Flows  measured  at  the  mouth  during  a  breath  are  read 
into  the  Ventilation  Model  and  the  output  from  the  Convection-Diffusion  Model,  the  exhaled 
gas  concentration  at  the  mouth,  can  be  compared  with  the  measured  concentrations. 


Figure  5-1.  Overview  of  the  Pulmonary  Model  Inputs  end  Outputs 
and  the  Validation  Process 
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The  Ventilation  Model  calculates  that  large  volume  and  ventilation  difierences  that  are 
created  by  gravity,  which  have  traditionally  been  called  “interr^onal  differences”  by 
physiol(^sts.'°  The  interregional  volume  and  ventilation  differences  are  a  result  of  die  stress 
field  in  the  lung  tissue  that  supports  the  weight  of  the  lung  tissue  and  blood.^  The 
Ventilation  Model  calculates  the  stress  and  strain  of  the  tissue  at  vertical  positions  during 
slow  and  rapid  breaths.  By  knowing  the  strain  of  the  tissue  at  different  lung  volumes,  the 
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regional  ventilation  pattern  can  be  determined.  Besides  gravity,  the  other  important 
parameters  in  the  Ventilation  Model  are  the  shape  of  the  lungs,  the  kinematics  of  the  pleural 
boundary,  the  blood  volume  and  distribution,  the  airway  lengths  and  diameters,  the  pressure 
losses  in  the  airways  and  the  elastic  characteristics  of  the  lung  parenchyma. 

The  Ventilation  Model  is  made  up  of  four  regions,  each  positioned  at  a  different  vertical 
location  as  shown  in  Figure  5-2.  Each  region  is  connected  to  the  trachea  by  a  symmetric 
system  of  conducting  airways.  The  position  of  the  first  three  airway  generations  is  shown  in 
Figure  5-2.  The  elastic  characteristics  of  the  tissue  in  each  region  are  developed  through  a 
tissue  model.  The  gas  pressures  in  each  region  are  found  by  calculating  the  pressure  losses 
in  the  gas  as  it  fiows  through  the  conducting  airways  of  that  region.  The  regional  flows  for 
a  given  flow  at  the  mouth  are  calculated  by  balancing  the  tissue  pressures,  the  gas  pressures 
in  the  airways  and  the  weight  of  the  tissue  and  blood  in  the  model.  The  output  from  the 
Ventilation  Model  is  the  four  regional  volumes  and  the  regional  flows  during  a  breath.  In 
order  to  validate  this  model  experiments  in  which  a  bolus  of  radioactive  xenon  gas  was 
inhaled  and  then  exhaled  have  been  simulated. 

The  Convection-Diffusion  Model  calculates  the  distribution  of  a  tracer  gas  along  the  axial 
pathways  in  the  lung  by  solving  the  one-dimensional  convection-diffusion  equation.  This 
model  is  also  comprised  of  four  regions  that  ventilate  in  parallel.  The  gas  volume  of  each  of 
the  conducting  airways  in  each  region  is  represented  physically  in  this  model  as  a  single 
trumpet  shaped  channel,  as  shown  in  Figure  5-3.  The  length  of  the  model  is  the  axial  length 
down  an  airway,  while  the  area  of  the  channel  is  the  total  cross-sectional  area  of  all  the 
airways  of  a  p^cular  airway  generation.  The  regional  volumes  and  bulk  flows  in  each 
region  during  a  breath  are  calculated  with  the  Ventilation  Model.  Using  these  flows  and 
volumes,  the  one  dimensional  convection-dififiision  equation  is  numerically  solved  in  each 
region  to  determine  the  specific  concentrations  of  a  single  tracer  gas  throughout  the  four 
regions  during  a  breath.  Single  breath  nitrogen  (Nj)  washout  tests  are  simulated  with  the 
model  in  order  to  validate  it. 

5.2  Ventilation  Model. 

5.2.1  Theory. 

The  lung  parenchyma  is  assumed  to  be  isotropic  and  homogeneous,  i.e.  the  tissue  and  gas 
form  a  continuous  material.  A  system  of  symmetric  conducting  airways  deliver  gas  to  the 
alveoli  in  each  r^on.  Normal  forces  at  the  pleural  surface  are  the  only  forces  supporting 
the  lung  weight  and  the  lung  is  free  to  slide  over  the  pleural  boundary.  In  order  to  simplify 
the  problem,  it  is  also  assumed  that  the  tissue  does  not  generate  a  shear  stress  to  oppose  a 
shear  strain. 
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Figure  5*2.  The  Four  Regiont  of  the  Pulmonary  Model. 
The  First  Three  Airway  Qenerations  (Gens  1  *3) 
and  the  Trachea  (Gen  0)  Are  Shown. 


A  free  body  diagram  of  a  thin  horizontal  section  of  the  idealized  lung  parenchyma  with  a 
thicimess  Az  is  show  in  Figure  S-4.  These  distance  from  the  lung  apex  to  top  surface  is  z,  A 
is  the  cross  sectional  area  of  the  section  at  z,  Ppl  is  the  pleural  pressure  at  z,  p  is  average 
density  of  the  homogeneous  material  in  the  section  and  g  is  the  acceleration  due  to  gravity. 

The  force  acting  on  the  upper  surface  is  A(P-Falv)  and  the  force  acting  on  the  lower  surface 
is 

A(/»-Palv)  +  — (i4(i»-f*a^))Az  .  [5-7] 

dz 


47 


Pleural  pressure  acts  on  the  outside  sur&ce  and  the  vertical  projection  of  this  fierce  is 


az 


[5*8] 


Excq>t  at  very  low  lung  volumes,  pleural  pressure  is  negative  and  this  force  supports  the 
lung  weight.  The  weight  of  the  section  is  pgAAZ.  The  inertia  of  the  tissue  and  the  viscous 
losses  that  occur  in  the  tissue  while  it  is  inflating  or  deflating  are  takra  to  be  n^gible. 
Therefore,  equilibrium  requires  that  the  forces  up  must  be  equal  the  forces  down.  Equating 
the  vertical  forces  in  Figure  S-4  and  simplifying  yields  the  following  expressimi 

—  [A(P-Palv)]  Az  +  iP-Palv)  ^  £iz  *  Ppl  ^  Az  *gpAAz=0  .  [5-9] 

dz  az  dz 

Static  equilibrium  at  the  pleural  surfoce  requires  that  Pp/=-(P-Rdv)  and  Equation  5-9 
becomes 


dz 


dPcdv 

dz 


=  -gp  • 


[5-10] 
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Figure  5-4.  Free-Body  Diagram  of  a  Thin  Horiiontal  Sliea  of  Lung  Tiaaua 


The  tissue  pressure  gradient,  dp/dz,  is  a  function  of  the  vertical  gradient  in  the  alveolar 
pressure  as  well  as  the  tissue  density.  If  Equation  5-10  is  int^rated  between  two  vertical 
positions  Zj  and  Zj,  then 

P^-  Pj^  f  pdz  *  (Paiv^  -  Patvp  ,  [5-11] 

where  the  intend  term  represents  the  pressure  difference  created  by  the  weight  of  die  tissue 
and  blood  wd  the  subscript  indicates  the  vertical  location. 

The  diffemice  in  the  alveolar  pressures  between  two  vertical  r^ons  is  obtained  by 
calculating  the  gas  pressures  in  the  bronchi  that  connect  the  two  r^ons.  Any  two  r^rnis  in 
the  lung  have  airway  that  meet  at  some  bronchus  which  is  called  the  meeting  bronchus.  As 
shown  in  Figure  5-5,  the  alveolar  pressure  of  each  region  minus  the  pressure  difference, 
APaw,  between  the  alveolar  gas  and  the  gas  pressure  at  the  me^g  bronchus  equals  the 
pressure  at  the  distal  end  of  the  meeting  bronchus 

Palv^~LPuw^•PalvJ~lPawJ  .  [5-12] 


In  Figure  5-5  the  Pbr  is  the  bronchial  pressure  and  the  n  subscript  is  the  generation  number. 
The  pressure  change  along  the  airway,  APaw,  equals  the  sum  of  the  pressure  losses  across 
each  bronchus  dP  in  the  pathway  b^ween  the  alveoli,  which  start  at  genmtion  17,  and  the 
meeting  bronchus 


49 


^PaWf  =  5^ 

»“16 


[5-13] 


where  the  pressure  change  across  each  bronchus  equals  the  pressure  at  the  distal  end  of  the 
bronchus  minus  the  pressure  at  the  proximal  end.  The  pressure  change  in  each  bronchus  is 
due  to  viscous  losses  (dPv)  and  convective  pressure  changes  {SPc)  as  the  gas  passes  from  one 
generation  of  airway  to  the  next,^‘ 

6Pbr  =  bPv*bPc  .  [5-14] 

When  Equation  S-12  is  substituted  into  Equation  S-11, 


P^  -  Pj  •  g  j  pdz  *  (AFow,  -  ^Pawp  •  [5-15] 

*1 

The  r^onal  volumes  and  flows  in  the  Ventilation  Model  are  calculated  by  satisfying 
Equation  5-15. 
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S.2.2  The  Lung  Model. 


As  shown  in  Figure  5-2,  each  lung  is  represented  as  four  regions  that  ventilate  in  parallel. 
Each  region  has  the  same  tissue  mass.  The  pressures  throughout  a  region  are  constant  except 
in  the  neighborhood  of  the  larger  bronchi  where  the  local  tissue  pressure  deprads  on  the 
relative  expansion  of  the  bronchi  compared  to  the  expansion  of  the  tissue  in  the  r^ion.  The 
alveolar  gas  volume  in  each  region  is  connected  to  the  trachea  by  a  homogeneous  system  of 
conducting  airways.  In  a  homogeneous  airway  system,  each  pathway  from  the  trachea  to  the 
alveolated  airways  is  identical.  The  gas  flow  through  each  pathway  within  a  region  is  taken 
to  be  the  same,  therefore,  the  pressure  change  APaw  is  the  same  for  all  airways  in  a  region 
and  the  alveolar  pressure  is  the  same  throughout  a  region.  Weibel’s  Symmetric  Lung  Model 
is  used  to  obtain  the  lengths  and  diameters  of  each  of  the  airway  generations  in  the 
Ventilation  Model.  In  Weibel’s  Symmetric  Lung  Model  A,  the  lung  is  composed  of  a 
symmetric  airway  system  of  23  generations.  Each  airway  division  distal  from  the  trachea  is 
labelled  as  a  generation.  Generation  1  is  the  two  airways  that  connect  the  right  and  left  lung 
to  the  trachea  and  generation  23  is  the  deepest  airway.  Alveoli  start  at  generation  17.  Table 
S-1  lists  the  lengths  and  diameters  of  the  airways  in  ftis  model. 

Figure  5-2  shows  the  lung  model  along  with  the  first  three  airway  generations.  Generations 

4- 16  are  all  located  in  the  region  where  the  third  generation  ends.  The  bronchial  pressures  in 
generations  1-16  are  calculated  in  the  Ventilation  Model.  The  gas  velocities  in  the  higher 
generations  are  too  low  to  produce  significant  pressure  changes.  All  four  regions  of  the 
model  meet  at  the  first  generation,  and  so  generation  1  is  the  meeting  bronchus  in  Equation 

5- 15. 


5.2.3  Parenchymal  Elastic  Recoil  Characteristics. 

A  model  of  the  lung  tissue  was  developed  as  part  of  the  Ventilation  Model  in  order  to  take 
into  account  the  elastic  characteristics  of  the  tissue,  the  shape  of  the  lung,  the  blood  volume 
in  the  lung,  the  general  elastic  characteristics  of  the  tissue  and  the  effect  of  airway  closure 
and  opening.  This  model  determines  the  tissue  elastic  characteristics  of  a  test  subject  from 
data  on  their  residual  volume  (RV),  total  lung  capacity  (TLC)  and  their  measured  lung 
pressure-volume  characteristics. 

The  basic  mechanical  unit  of  the  tissue  model  consists  of  a  single  terminal  bronchus 
(generation  17)  which  services  a  respiratory  unit  (generations  18-23).  The  mechanical  unit 
represents  the  gas  exchanging  portion  of  the  lungs  which  comprises  more  than  90%  of  the 
gas  volume  of  the  lung.  All  respiratory  units  in  the  tissue  model  are  identical.  The 
respiratory  unit  can  remain  patent  until  zero  gas  volume,  but  the  terminal  bronchus,  which 
is  the  sight  of  airway  closure  in  the  lungs,^^  acts  as  a  valve  which  is  either  fiilly  open  or 
closed.  A  pressure-fractional  volume  curve  for  the  respiratory  unit  (P-FVru)  could  be 
constructed  by  expanding  the  unit  from  the  completely  collapsed  state,  zero  gas  volume,  to 
its  maximum  volume  and  then  deflating  the  unit  back  to  zero  volume.  The  factional  volume 
(FV)  is  the  volume  of  a  lung  unit  normalized  with  respect  to  its  maximal  volume.  The  gas 


Table  5-1.  Weibel’s  Symmetric  Lung  Model  A, 
Total  Volume  of  Lung  Model  (Vm)  is  4820  cc. 


Airway 

Number  of  Airways 

Generation 

Generation 

Alveolar  V<^nme 

Oeoeralioo 

per  Oeuention 

Length 

Diameter 

per  Geoeratioo 

(NPO) 

1  (cm) 

(on) 

Valv(oc) 

16 


32 


64 


128 


256 


512 


1024 


2048 


4096 


8192 


16384 


32768 


65536 


131072 


262144 


524288 


1048576 


2097152 


4194304 


8388608 


4.76 


1.90 


0.76 


1.27 


1.07 


0.90 


0.76 


0.64 


0.54 


0.46 


0.39 


0.33 


0.27 


0.23 


0.20 


0.165 


0.141 


0.117 


0.099 


0.083 


0.070 


0.059 


0.050 


1.22 


0.83 


0.56 


0.45 


0.35 


0.28 


0.23 


0.186 


0.154 


0.130 


0.109 


0.095 


0.082 


0.074 


0.054 


0.050 


0.C47 


0.045 


0.043 


0.041 


0.041 


6. 


22.2 


63.2 


222.2 


444.5 


892.2 


1523.0 
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pressure  measured  at  a  volume  when  the  lung  unit  is  sealed  reflects  the  tension  in  the  tissue 
and  is  called  the  elastic  recoil  pressure  (P)  or  simply  the  unit  pressure.  The  P-FVru  curves 
used  in  the  tissue  model  are  obtained  by  extrapolating  the  upper  half  of  a  lung  pressure- 
volume  curve  (P-FVlung)  to  zero  gas  volume.  At  a  fractional  gas  volume  greater  than  0.50 
all  of  the  airways  are  presumed  open  and  the  lung  pressure-fractional  volume  curve  and 
tissue  fractional  volume  curve  should  be  similar.  A  cubic  curve  is  used  to  describe  the  upper 
half  of  the  lung  pressure-fractional  volume  curve  and  a  straight  line  is  used  to  describe  the 
lower  half,  as  given  in  the  following  equation. 


FVnuiiF)  =  a(P-Po.o)  ' 

P  =  0  for  OsFKnitt^O.5 
P  =  16  a  —  for  0.5<FV'iii«sl.0 


The  value  of  a  for  a  subject  is  determined  from  the  linear  part  of  his  P-FVlung  curve,  which 
is  between  approximately  20%  and  80%  of  the  vital  capacity  (VC).  The  value  of  Po.s  equals 
the  transpulmonary  pressure  at  FVlung  =  0.5  during  inhalation  and  exhalation.  The  dashed 
line  in  Figure  5-6  shows  the  P-FVru  relation  for  a  subject  with  P05  =  5  cmHjO  during 
exhalation,  a  =  0.05  cmH20  and  a  hysteresis  of  4  cmH20  during  a  VC  breath.  Both  the 
inhalation  and  exhalation  curves  are  taken  to  have  the  same  shape,  only  they  are  separated  by 
the  hysteresis. 

The  state  of  each  terminal  bronchus,  whether  it  is  open  or  closed,  depends  on  the  transmural 
(Ptm)  that  acts  on  it.  The  transmural  pressure  is  defined  as  the  intrabronchial  gas  pressure 
(Pbr)  minus  extrabronchial  pressure  (Px),  the  pressure  acting  on  the  outside  wall.  Since  the 
terminal  bronchi  are  deeply  embedded  in  the  tissue,  the  extrabronchial  pressure  equals 
pleural  pressure  under  the  no  fiow  condition,  so  the  transmural  pressure  acting  on  the 
terminal  bronchi  equals  pleural  pressure  under  no  flow  conditions. 

The  transmural  pressure  at  which  a  terminal  bronchus  closes  is  called  the  critical  closing 
pressure  (Pcc).  The  amount  of  gas  trapped  and  the  minimum  gas  volume  of  each  respiratory 
unit  is  determined  by  the  critical  closing  pressure  of  its  terminal  bronchus.  Once  closed  each 
terminal  bronchus  requires  a  significantly  higher  transmural  pressure  before  reopening.  The 
transmural  pressure  at  which  the  terminal  bronchus  reopens  is  called  the  critical  opening 
pressure  (Pco).  Unfortunately,  there  is  little  data  available  on  the  distribution  of  the  critical 
pressures  so  it  is  taken  to  have  a  normal  distribution.  The  mean  critical  opening  and  closing 
pressure  of  the  probability  density  function  are  Pcc  and  Pco  and  the  stand^d  deviations  of 

the  distributions  are  (7cc  and  (Tco.  The  sequence  of  the  terminal  bronchi  opening  is  taken  to 
be  the  opposite  of  the  closing  sequence,  i.e.  the  first  airways  to  close  are  the  last  to  open. 

The  maximum  range  of  the  critical  pressures  has  been  estimated  to  be  4.0  cmH20^^  and  the 
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Figure  5-6.  The  Pressure-Fraction  Volume  Relation  of  a  Mechanical  Unit  (P-FVmu) 
Is  Shown  as  a  Solid  Line  and  for  a  Respiratory  Unit  (P-FVru)  as  a  Dashed  Line 
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standard  deviations  of  the  critical  pressures  are  (Tcc  =  0,5  and  (Jco  =  0,8.  Figure  5-7 
shows  examples  of  the  probability  density  functions  for  the  critical  opening  and  closing 
pressures. 

Figure  5-7.  Example  of  Probability  Density  Function  for  Critical  Closing/Opening 
Pressures  of  Terminal  Bronchi  Used  in  Tissue  Model 
Pcc  =  -0.8  cmHgO.  Pco  =  4.2  cmHjO 


-2  -1  01234567 

Critical  Transmural  Pressure  (cm  H2O) 


The  pressure-fractional  volume  relation  of  the  mechanical  unit  (P-FVmu)  represents  the 
elastic  characteristics  of  the  respiratory  unit  and  its  terminal  bronchus.  The  P-FVmu  curve  is 
identical  to  the  P-FVru  curve  as  long  as  the  terminal  bronchus  remains  open.  Once  the 
terminal  bronchus  closes,  the  mechanical  unit  is  incompressible  and  any  further  lowering  of 
the  tissue  pressure  causes  no  decrease  in  gas  volume.  The  P-FVmu  of  a  mechanical  unit 
during  a  vital  capacity  breath  is  shown  in  Figure  5-6  as  a  solid  line.  The  terminal  bronchus 
of  this  mechanical  unit  has  a  Pcc  =  0.0  cmHjO  and  Pco  =  5  cmHjO.  During  exhalation, 
once  the  terminal  bronchus  closes,  the  mechanical  unit  is  incompressible  and  the  fractional 
volume  remains  at  0,2.  When  the  tissue  pressure  is  increased  to  5  cmHjO,  the  terminal 
bronchus  reopens  and  the  mechanical  unit  rapidly  expands  to  its  new  equilibrium  state  on  the 
P-FVru  curve. 
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The  fractional  gas  volume  of  the  tissue  at  a  pressure  P  equals  the  sum  of  the  fractional  gas 
volumes  of  all  the  mechanical  units  that  comprise  the  tissue.  The  tissue  pressure-fractional 
volume  (P-FVtis)  relations  are  obtained  by  accounting  for  the  closing  and  subsequent  opening 
of  the  terminal  bronchi  in  all  the  mechanical  units  of  the  tissue.  The  P-FVtis  curve  that 
results  when  the  P-FVru  curve  in  Figure  S-6  is  combined  with  the  opening  and  closing 
distribution  in  Figure  5-7  is  shown  in  Figure  5-8.  The  code  used  to  calculate  the  P-FVtis 
curve  is  given  in  Appendix  C-1. 

The  total  tissue  mass  of  the  lungs  is  taken  to  be  a  function  of  a  subject’s  TLC.  At  TLC  all 
regions  of  the  lung  are  fully  expanded  and  have  the  same  density  of  0.12  gm<cm'^  per  cc. 

RV  varies  from  subject  to  subject  and  in  the  same  subject.  For  each  subject  the  ratio  of  RV 
divided  by  TLC  and  the  tissue  minimum  fractional  volume  (FVti^niJ  are  important 
parameters  in  calculating  the  lung  pressure-fractional  volume  relation  at  low  lung  volumes 
because  these  parameters  determine  the  fraction  of  the  respiratory  units  that  are  closed  at  low 
lung  volumes.  In  order  to  quantify  the  relation  between  RV/TLC  and  FVtiSi^,  a  new 
variable  is  introduced,  the  difference  between  the  minimal  fractional  volume  of  the  lung  and 
the  tissue, 


t  =  (^)  -  FVtis^  [5-17] 

The  graphical  meaning  of  \l/  is  shown  in  Figure  5-8,  where  the  P-FVtis  and  the  calculated 
lung  Pressure-Frictional  Volume  (P-FVlung)  curve  for  a  subject  are  shown.  For  this  subject, 
RV/TLC  =  0.26  and  FVtis„un  =  0-21,  making  ^  =  0.05.  If  all  of  the  lung  tissue  at  RV  is 
not  at  the  minimal  fractional  volume  then  \p>0,  which  is  the  case  for  most  young,  healthy 
subjects.  If  all  of  the  lung  tissue  is  at  the  minimal  fractional  volume  then  ^  =  0.  Thus 
can  only  be  equal  to  or  greater  than  zero. 

There  are  no  experimental  estimates  of  ip,  but  measurements  of  the  fraction  of  airways  closed 
at  RV  as  a  function  of  vertical  position  have  been  made.^^  By  choosing  \p  equal  to  0.05,  the 
fraction  of  airways  closed  along  the  vertical  axis  of  the  lung  model  was  similar  to  that 
measured  in  these  experiments. 

The  P-FVlung  relation  and  the  degree  of  airway  closure  are  determined  by  creating  a  lung 
with  a  physical  geometry  like  the  human  lung,  a  ellipsoidal  shape,  and  solving  Equation  5-10 
with  the  dveolar  pressure  gradient  set  equal  to  zero. 


dz 


=  -gP 


[5-18] 


the  hydrostatic  pressure  equation.  By  filling  the  ellipsoidal  shaped  lung  model  with  the 
imaginary  tissue  and  solving  the  hydrostatic  pressure  equation  Ae  P-FVlung  can  be 
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Figure  5-8.  The  Relation  Between  tp  and  the  P-FVtis  and  P-FVIung  Curves. 
P-FVtis  Is  Shown  as  a  Dashed  Line  and  P-FVIung  is  Shown  as  a  Solid  Line. 
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calculated.  The  P-FVlung  relationships  are  obtained  by  dividing  the  ellipsoidal  sluq)ed  model 
into  fifty  thin  sections  and  numerically  solving  the  hydrostatic  pressure  equation  at  the 
bounda^  of  each  slice.  The  lung  model  then  performs  a  VC  inhalation  and  exhalation,  while 
the  elastic  recoil  pressure  is  measured  at  a  fixed  vertical  position  that  is  similar  to  the 
position  of  an  esophageal,  balloon.  In  doing  this  calculation  the  weight  of  the  blood  is  also 
included  in  the  density  term.  The  code  for  calculating  the  P-FVlung  relation  is  given  in 
Appendix  C-2. 

5.2.4  Blood  Volume  Distribution 

The  blood  volume  and  distribution  affects  the  tissue  pressure  since  supporting  stresses  in  the 
tissue,  the  dp/dz  term  in  Equations  S-10  and  5-15,  must  support  the  weight  of  the  blood  as 
well  as  the  tissue  weight.  Information  on  the  volume  distribution  of  the  blood  in  the 
pulmonary  capillaries  is  difficult  to  obtain,  and  to  the  best  of  our  knowledge,  no  data  is 
available  on  humans,  although  there  have  been  estimates^  of  total  blood  volume  ranging 
from  approximately  100  cm^  to  approximately  200  cm^  In  order  to  calculate  the  blood 
volume  distribution  along  the  vertical  axis,  the  pulmonary  capillaries  are  viewed  as  sheets 
that  surround  the  alveoli.^  The  sheet  area  per  gram  of  tissue  is  taken  as  a  constant 
throughout  the  lung  and  is  estimated  to  be  100  square  centimeters  per  gram  of  tissue.  The 
height  of  the  capillaries  (h)  is  calculated  using  experimental  results.^^’^  In  these  studies  it 
was  found  that  the  capillary  height  in  cats  is  a  linear  function  of  the  capillary  tiansmural 
pressure,  the  pressure  differences  between  the  capillary  pressure  (Pcap),  and  the  alveolar  gas 
pressure. 


h  =  +  h'(Pcap  -  Pah) 


[5-19] 


where 


\  =  4.3  X  lO-^cw 


[5-20] 


and 


h'  =  22  X  lO  ’cm  H^O  .  [5-21] 

This  relationship  holds  for  Pcap  greater  than  Palv,  but  for  negative  transmural  pressures  the 
capillary  is  collapsed  and  h  =  0.  The  capillary  blood  pressure  is  the  sum  of  the  arterial 
pressure  generate  by  the  heart  and  the  hydrostatic  pressure  that  results  in  the  change  in 
vertical  position  as  the  blood  moves  from  the  heart  through  the  lungs.  Figure  5-9  shows  the 
vertical  distribution  in  one  of  the  lungs  of  the  tissue  model  at  60%  VC.  The  model  predicts 
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that  the  blood  volume  will  increase  as  lung  volume  increases,  an  increase  caused  by  the 
downward  movement  of  the  base  of  the  lung  relative  to  the  heart. 


Figure  5-9.  The  Summed  Blood  Volume  Moving  From  the  Lung  Apex  to  the  Base 
in  a  Single  Lung.  Total  Pulmonary  Blood  Volume  Would  Be  Twice  This  Amount. 
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The  tissue  model  is  also  used  to  obtain  the  Pressure  Fractional  Volume  relations  for  each  of 
the  regions  in  the  Ventilation  Model  (P-FVreg).  These  are  obtained  by  dividing  the  fifty 
horizontal  sections  of  the  tissue  model  into  four  regions  of  equal  mass.  The  center  of  gravity 
(CG)  of  each  of  these  four  regions  is  determined  and  the  tissue  model  is  taken  through  a  VC 
inhalation  and  exhalation  while  the  elastic  recoil  pressure  at  the  regional  CG’s  and  the  gas 
volume  of  the  r^ion  are  recorded.  The  other  variable  measured  is  the  difference  in  the 
elastic  recoil  pressure  between  the  regional  CG’s,  which  takes  into  account  the  weight  of  the 
tissue  and  blood  between  the  regional  CG’s.  The  P-FVreg  curves  represent  the  elastic 
characteristics  of  each  of  the  four  regions  in  the  Ventilation  Model. 

5.2.5  Transmural  Pressures  Acting  on  the  Bronchi 

The  bronchi  in  the  Ventilation  Model  are  treated  as  compliant  tubes  whose  diameters  are  a 
function  of  the  transmural  pressure  acting  on  their  walls.  The  transmural  pressure  (Ptm) 
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equals  the  gas  pressure  inside  the  bronchus  (Pbr)  minus  the  extrabronchial  pressure  (Px),  the 
pressure  acting  on  the  outside  wall  of  the  bronchus. 


Ptm  =  Pbr  -  Px 


[5-22] 


The  tissue  pressure  throughout  each  region  of  the  model  is  equal  to  P  excq>t  in  the 
neighborhood  of  the  larger  bronchi.  The  extrabronchial  pressure  is  the  sum  of  the  local 
tissue  pressure  and  the  alveolar  pressure.  Px  depends  on  how  the  brcmchus  is  expanded 
relative  to  the  tissue  in  the  surrounding  lung  r^on  since  there  is  a  mechanical 
interdependence  between  the  bronchus  and  the  tissue  that  surrounds  it.^’”  If  the  bronchus 
expands  uniformly  with  the  surrounding  region,  then  the  extrabronchial  pressure  equals  the 
pleural  pressure.  In  this  case,  the  tissue  around  the  bronchus  expands  as  if  the  space 
occupied  by  the  bronchus  was  filled  with  tissue. 

Whoi  a  bronchus  does  not  expand  uniformly  with  the  surrounding  tissue,  the  local  tissue 
pressure  does  not  equal  the  r^<Hial  tissue  pressure  and  Px  no  longer  equals  pleural  pressure. 
For  example,  during  a  rsq)id  exhalation,  the  bronchi  defiate  faster  than  the  surrounding  tissue 
because  they  can  be  exposed  to  large  n^ative  Ptm’s.  When  this  hjqrpois,  Px  becomes 
greater  than  the  pleural  pressure  because  of  the  mechanical  interdepradence  and  this  increase 
in  Px  aids  in  keying  the  bronchi  from  collapsing.  The  effect  of  bronchial  and  tissue 
interdependoice  in  the  Ventilation  Model  is  assumed  to  be  important  during  exhalation,  it  is 
taken  to  be  insignificant  during  inhalation.  The  magnitude  of  the  mechanical  interdqpendence 
has  been  estimated  using  the  interdependence  parameter  that  was  introduced  to  explain  some 
experimental  data.^'  The  transmural  pressure  is  calculated  by  solving  the  following  equation. 


Ptm  ^  P  -  (^APow 


[5-23] 


where  P  is  the  regional  elastic  recoil  pressure,  APaw  is  the  total  pressure  change  betweoi 
the  alveoli  and  that  airway,  and  K  is  the  interdependence  parameter  which  ranges  from  one 
to  approximately  three. 

The  elastic  properties  of  the  bronchi  and  the  tissue  have  been  shown  to  be  similar.” 
Therefore,  the  bnmchial  diameters  and  lengths  are  treated  as  a  functirm  of  the  gas  volume  of 
the  lung  V  and  varied  with  the  cube  root  of  the  fractional  gas  volume  of  the  lung 
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where  I  is  the  length  and  d  is  the  diameter  of  an  airway  generatitm  and  and  are  the 
maximum  airway  lengths  for  a  generation  at  TLC.  Equation  5-2A  is  used  to  solve  for  the 
relative  lengths  and  diameters  of  the  bronchi. 

During  exhalation  the  bronchi  are  exposed  to  negative  transmural  pressures.  Unfortunately, 
there  are  limited  data  on  the  diameters  changes  of  the  bronchi  under  negative  transmural 
pressures.  During  exhalation  the  conducting  airway  system  and  latex  tubes  both  exhibit  flow 
limitation,  that  is  the  flow  approaches  an  asymptote  (for  a  constant  lung  volume)  as  the 
applied  pressure  is  increased.  The  relation  of  the  airway  diameter  d  to  a  negative  transmural 
pressure  is  described  by  the  modified  tube  law,^ 


d^^d^iX-^)  [5-25] 

where  do  is  the  airway  diameter  of  a  generation  at  Ptm  =  0  and  the  coefficient  e  is 
proportional  to  the  bending  stiffness  of  the  bronchial  wall  at  Ptm  =  0. 

5.2.6  Solution  Procedure. 

The  regional  fiows  are  calculated  for  a  given  flow  at  the  mouth  Qm  at  a  given  lung  volume. 
The  calculation  is  made  simultaneously  solving  Equation  5-15  and  equation  of  conservation 
of  mass, 

=  +  [5-26] 

where  the  Qi,  Q2,  Qa,  and  Q4  are  the  regional  fiows.  The  calculation  of  the  regional  flows  is 
done  numerically.  The  breath  is  divided  into  small  time  increments  of  size  At  and  the 
regional  flows  are  calculated  at  the  end  of  each  time  increment  by  determining  the  regional 
flows  that  satisfy  Equation  5-15  and  the  conservation  of  mass,  Equation  5-26.  The  c^e  for 
the  Ventilation  Model  is  given  in  Appendix  C-3. 

5.2.7  Model  Results  and  Validation. 

Figure  5-10  shows  the  volumes  of  the  four  regions  of  the  Ventilation  Model  during  a  slow 
VC  inhalation  by  a  subject  with  a  RV  of  2  liters  and  a  TLC  of  7.5  liters.  The  tissue  P-FVtis 
and  lung  P-FVlung  relation  for  this  subject  are  shown  in  Figure  5-8.  The  distribution  of 
opening  and  closing  pressures  for  this  subject  are  shown  in  Figure  5-7. 

Figure  5-10  provides  a  good  overview  of  the  interregional  ventilation  differences  during  a 
slow  breath.  The  VC  for  this  subject  was  5.6  liters.  Most  of  the  gas  inhaled  near  RV  enters 
the  apical  regions.  Regions  1  and  2,  and  the  volume  curves  for  these  two  regions  have  the 
steepest  slope.  Region  4  undeigoes  the  smallest  expansion  during  the  inhalation  of  the  first 
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liter  of  gas  because  it  contains  such  a  large  fraction  of  closed  terminal  bronchi.  Once  all  the 
terminal  bronchi  open,  the  basilar  regions  become  more  compliant  and  take  in  a  larger 
fraction  of  the  inhaled  gas,  as  indicated  by  the  steeper  slope  of  the  r^ional  volume  curves 
for  Regions  3  and  4  after  1.5  liters  of  gas  have  been  inhaled.  This  causes  the  regional 
volume  differences  to  decrease  as  lung  volume  increases.  At  TLC,  the  four  regions  contain 
almost  the  same  amounts  of  gas. 

The  emptying  of  the  four  regions  is  also  sequential,  the  slope  of  the  volume  curve  for 
Regions  1  and  2  becomes  progressively  steeper  as  the  gas  volume  decreases,  the  two  apical 
regions  contribute  more  gas  to  the  expiate  as  lung  volume  decreases.  Initially  the  basilar 
regions  contribute  the  largest  fraction  of  gas  and  their  gas  volumes  decrease  the  fastest.  As 
lung  volume  approaches  RV,  the  terminal  bronchi  in  the  basilar  regions  begin  to  close, 
reducing  the  flow  from  these  regions  and  the  regional  volume  curve  for  Region  4  becomes 
flatter.  At  the  end  of  exhalation,  almost  all  of  the  decrease  in  total  gas  volume  occurs  in  the 
upper  two  regions  because  of  the  laige  fraction  of  closed  terminal  bronchi  in  the  basilar 
regions.  This  plot  of  the  volume  change  of  regions  at  different  vertical  positions  agrees  with 
experimentally  measured  curves  of  regional  volume  changes.*’ 

The  Ventilation  Model  has  been  used  to  calculate  both  the  inhaled  distribution  and  the 
subsequent  washout  of  radioactive  xenon  (Xe)  gas.*^  *^  In  these  experiments,  a  subject 
inhales  a  bolus  of  radioactive  Xe  gas  at  different  lung  volumes  and  at  different  inh^tion 
flow  rates.  At  the  end  of  inhalation,  the  subject  performs  a  breath  hold  and  the  vertical 
concentration  of  Xe  gas  in  the  lungs  is  measured.  During  the  subsequent  exhalation,  the  Xe 
concentration  at  the  mouth  is  measured. 

In  order  to  perform  an  accurate  simulation  of  the  Xe  washin  experiments,  the  transport  of  the 
bolus  through  the  upper  airways,  the  mouth,  the  pharynx  and  trachea  is  analytically  modeled. 
This  is  done  by  representing  the  upper  airways  as  a  pipe  with  a  cross-sectional  area  of  2.5 
cm^  and  a  length  of  25  cm.  A  Xe  bolus  of  5  cm^  is  placed  at  the  pipe  entrance  and 
convected  down  the  pipe  with  a  uniform  velocity.  The  mixing  in  the  pipe  is  calculated  by 
using  an  effective  dispersion  coeflScient.  The  bolus  is  mixed  with  air  as  it  moves  through  the 
upper  airways  and  by  the  time  it  reaches  the  end  of  the  pipe,  approximately  125  cm^  of  gas 
contains  radiation  in  the  slow  inhalation.  The  end  inhalation  Xe  concentrations  calculated 
with  the  Ventilation  Model  reflect  how  the  model  expands  during  the  inhalation  of  this  125 
cm^  of  gas.  During  the  subsequent  exhalation,  the  washout  of  the  Xe  at  the  mouth  is 
recorded  as  a  function  of  the  exhaled  volume.  If  there  are  Xe  concentration  differences  in 
the  lung  at  the  end  of  inhalation,  the  Xe  concentration  at  the  mouth  reflects  the  emptying 
pattern  of  the  different  lung  regions  during  exhalation. 

The  results  of  the  inhalation  simulations  are  shown  in  Figures  5-11  and  5-12  for  slow 
inhalations  and  fast  inhalations.  The  results  are  expressed  as  Xe  concentration  versus  z,  the 
vertical  height  of  the  lung.  The  concentrations  are  shown  for  an  Xe  bolus  inhaled  at  RV  in 
Figure  5-11  and  at  55%  VC  in  Figure  5-12.  The  Xe  that  enters  each  region  is  assumed  to  be 
completely  mixed  in  the  alveolar  gas  of  the  region.  The  Xe  concentration  curves  for  each 
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region  of  the  model  are  horizontal  lines  when  plotted  against  z,  thus  the  step  sh£q)ed 
appearance  of  the  curves  in  Figures  5-11  and  5-12. 

The  agreement  between  the  calculated  and  measured  end  inhalation  Xe  concentrations  is 
good.  Both  the  measured  and  calculated  results  show  that  as  the  inhalation  flow  rate 
increases,  the  distribution  of  the  Xe  along  the  vertical  axis  becomes  more  uniform.  In  the 
Ventilation  Model,  the  fraction  of  the  flow  going  to  the  upper  half  near  RV  decreases  as  the 
inhalation  flow  rate  increases.  At  55  %  VC  the  measured  and  calculated  end  inhalation 
concentrations  are  also  similar.  Both  the  measured  and  calculated  Xe  concentrations  show 
that  as  the  flow  rate  is  increased,  the  distribution  of  the  inhaled  bolus  in  the  lung  becomes 
more  uniform.  The  good  agreement  between  the  calculated  and  measured  results  indicates 
that  the  Ventilation  Model  accurately  describes  the  filling  pattern  of  the  lung  near  RV  and 
55%  VC  for  different  flow  rates. 

The  measured  single  breath  washout  curve  for  a  rapid  exhalation  starting  at  TLC  were  not 
available  in  the  study  that  investigated  the  inhalation  distribution  of  the  Xe  bolus.  A  similar 
study  did  publish  Xe  single  breath  washout  curves  for  different  exhalation  flows,  but  no  data 
was  given  on  the  lung  volumes  of  the  subjects  or  the  exact  flows  during  the  washout.*’ 
Therefore,  slow  and  rapid  washout  tests  have  been  simulated  using  the  lung  volume  and  the 
regional  pressure-fractional  volume  relations  of  the  subject  whose  lung  volumes  were  used  to 
produce  the  simulations  shown  in  Figures  5-11  and  5-12.  The  calculated  Xe  washout  curves 
are  shown  in  the  top  part  of  Figure  5-13.  The  slow  washout  has  a  flow  of  ^^roximately 
200  cm’  •  sec  *  and  the  rapid  washout  has  a  flow  of  approximately  2000  cm’  •  sec  *.  This  flow 
is  reduced  as  a  lung  volume  approaches  RV  in  order  to  keep  the  pressure  changes  in  the 
conducting  airways  in  the  physiological  range.  The  Xe  concentrations  in  the  first  750  cm’  of 
exhaled  gas  are  not  shown  since  the  model  only  represents  the  emptying  of  the  well  mixed 
alveolar  gas,  not  the  initial  gas  where  difiusion  plays  a  role.  The  regional  Xe  concentrations 
at  TLC  are  set  up  in  the  simulations  and  experiments  by  slowly  inhaling  a  bolus  at  RV, 
which  places  most  of  the  Xe  in  the  apical  regions  (see  Figure  5-11).  The  washouts  are 
expressed  as  a  function  of  exhaled  volume,  which  is  expressed  as  a  %  VC. 

During  the  slow  washout,  the  calculated  Xe  concentration  at  the  mouth  slowly  increases,  thra 
shows  a  rapid  rise  once  approximately  80%  of  the  VC  has  been  exhaled.  This  rise  at  the  end 
of  the  washout  is  caused  by  the  closure  of  airways  in  the  basilar  r^ons  and  the  resulting 
decrease  in  compliance.  As  the  compliance  of  the  lower  regions  decreases  the  upper  regions 
provide  more  of  the  gas  to  the  expirate  which  causes  the  concentration  of  Xe  at  foe  mouth  to 
increase.  Whoi  foe  exhalation  rate  is  rapid,  foe  model  predicts  foe  lung  will  empty  more 
uniformly,  foe  Xe  concentration  at  foe  mouth  is  more  uniform  and  actually  takes  a  dip  at 
approximately  75%  VC.  At  foe  very  end  of  the  rapid  exhalation,  a  rise  in  foe  Xe 
concentration  occurs. 
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Figure  5-13.  Top:  Calculated  Xe  Washout  Curves  During  Slow/Fast  VC  Exhalation. 
Bottom:  Measured  Xe  Washout  Curves  During  Slow/Fast  VC  Exhalation. 
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The  measured  curves  are  shown  in  the  lower  part  of  Figure  5-13.  The  relative  positions  of 
the  measured  slow  and  rapid  wash  out  curves  are  similar  to  the  relative  positions  of  the 
calculated  slow  and  rapid  wash  out  curves.  The  emptying  pattern  of  the  Ventilation  Mottel 
appears  to  qualitatively  represent  the  emptying  pattern  of  human  lungs  during  the  rapid 
exhalations. 

5.3  The  Convection-Diflhsion  Model. 

5.3.1  Geometry  of  the  Model. 

The  dimensions  of  the  trumpet  channel  are  obtained  from  the  Weibel’s  Symmetric  Lung 
Model  A.  The  anatomical  data  is  given  as  airway  length  per  generation  (4)>  diameter  of 
each  generation  (dj,  number  of  airways  per  generation  (nawj  and  the  alveolar  volume  per 
generation  (Valvj  (see  Table  5-1).  The  last  2.2  cm  of  the  trumpet  channel  that  rq}resents 
Weibel’s  Symmetric  Lung  Model  A  is  shown  in  Figure  5-14.  The  Y  axis  is  the  tc^  cross- 
sectional  area  of  all  the  airways  and  the  X  axis  shows  the  distance  from  the  trachea.  The 
most  important  aspect  of  the  trumpet  channel  in  Figure  5-14  is  the  large  cross-sectional  area 
at  the  distal  end  of  the  trumpet.  Over  90%  of  the  volume  is  contained  in  the  last  centimeter 
of  the  channel.  The  cross-hatched  area  indicates  the  alveolar  volume  in  each  goieration. 

The  alveoli  are  not  situated  on  the  conducting  airways  as  individual  sacs,  instead  they  are 
found  in  groups  with  each  group  being  connected  to  the  conducting  airway  by  a  structure 
called  a  sac.^  Figure  5-15  shows  a  schematic  drawing  of  a  sac  with  its  attached  alveoli, 
which  is  called  an  alveolar  unit.  All  of  the  alveoli  in  the  transport  modd  are  assumed  to  be 
arranged  in  alveolar  units.  There  area  of  the  sac  opening  is  taken  to  equal  0.Q52  mmP,  the 
depth  of  the  sac  and  its  attached  alveoli  is  0.32  mm,  the  volume  of  the  sac  is  0.00454  mm’ 
and  the  volume  of  each  alveolus  is  0.00392  mm’.^  If  all  of  the  alveoli  are  assumed 
connected  to  sacs  there  are  4.65  alveoli  per  sac  and  the  volume  of  the  alveolar  unit  equals 
0.02671  mm’.  The  total  area  of  the  sac  openings  per  generation  (Aso)  equals  the  sac 
opening  area  multiplied  by  the  number  of  sacs  per  generation. 

The  anatomical  lung  data  in  Wdbel’s  Symmetric  Lung  Model  A  represoits  a  lung  volume  of 
4,820  cm’  and  this  data  is  scaled  to  the  starting  lung  volume  of  the  subject  whose  breath  is 
bdng  simulated.  The  scaling  is  performed  assuming  the  lengths  changes  are  proportional  to 
the  cube  toot  of  the  volume  change.”  The  reference  dimensions  define  the  anatomical  state 
of  each  regicm  at  the  start  of  inhalation.  During  a  breath  each  region  infiates  and  defiates 
uniformly  in  three  dimensions.  The  breathing  process  can  be  pictured  as  the  movemrat  of 
air  into  and  out  of  each  r^onal  channel  as  the  charmel  shown  in  Figure  5-3  expands  and 
contracts.  The  dimensions  of  the  channel  at  each  new  regional  volume  are  obtained  by 
multiplying  the  reference  dimensions  by  the  correct  power  of  the  cubic  scaling  factor. 


Figure  5-14.  Distance  From  the  Trachea  vs.  the  Log  of  the  Total  Cross-Sectional 
Area  Over  the  Distal  Most  2.2  cm  of  the  Model. 
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where  Vi  is  the  volume  of  the  ith  region  and  is  the  reference  volume  of  the  ith  region. 
Therefore,  at  a  given  volume  Vi  with  a  scaling  factor  s,  the  length  per  generation  is  sl^,  the 
diameter  is  sd„,  the  area  of  the  channel  is  s^A,,  the  alveolar  volume  per  generation  is  s^Valv^ 
and  the  area  of  the  sac  openings  per  generation  is  s^Aso,,  where  the  n  subscript  indicates  the 
airway  generation  number. 

Figure  5-1 5.  An  Alveolar  Unit  Consists  of  a  Sac  and  Its  Attached  Alveoli. 

The  Sac  Connects  the  Alveoli  to  the  Conducting  Airway. 


Alveolar  Unit 


Conducting  Airway 


5.3.2  Convection-Difiiision  Transport  Equation. 

The  mass  balance  for  a  tracer  gas  in  a  channel  with  a  cross-sectional  area  of  s^A  can  be 
established  by  considering  a  differential  element  of  the  trumpet  channel  with  the  length  Ax  as 
shown  in  Figure  5-16. 

The  rate  of  change  of  tracer  gas  in  the  element  is  given  by 

-^(sMCAx)  ,  [5-28] 

dt 

which  must  equal  the  net  gas  flux  of  tracer  gas  into  the  element.  Mathematically  this 
corresponds  to  the  following  equation 


4(5M(ac)C(jc)Aar)  = 
dt 


-^s^i4(*)/(jc)  Aac)  -s^Aso(fi)Jdlvix)  An 
dx 


[5-29] 


where  t  represents  time,  C  the  airway  concentration,  J  the  axial  gas  flux  in  the  airway,  Jalv 
is  the  gas  flux  into  the  alveolar  unit,  n  is  the  generation  number  and  An  equals  Ax  divided  by 
s/,  the  fraction  of  an  airway  generation  in  the  Ax  segment.  For  convenience,  time  has  not 
been  listed  as  an  independent  variable  in  Equation  5-29,  but  it  should  be  noted  that  C,  J, 

Jalv,  and  s  are  all  functions  of  time.  The  flux  in  the  channel  is  given  by 

J(x)  =  [/(x)C(x)  [5-30] 

dx 

where  the  axial  velocity 

U(x)  .  ,  [5-31] 

sMcx 
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Qalv  = 


dValv 

dx 


[5-36] 


To  make  the  problem  tractable,  two  assumptions  are  made  concerning  the  tracer  gas  flux  into 
the  alveolar  unit.  First,  the  radial  dependence  is  eliminated  by  approximating  the  radial 
concentration  gradient  at  the  sac  opening  with 
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[5-37] 


Cah-C 

4(|+SD)/2 


Secondly,  the  concentration  at  the  sac  opening  is  assumed  to  be  the  average  of  the  airway 
concentration  C  and  the  alveolar  concentration  Calv,  Cso  =  (Calv+C)/2.  Substituting  these 
equations  into  Equation  5-16,  the  transport  into  the  alveolar  region  is 


s^AsoJalv  =  —is^ValvCalv) 
dt 


ggl„(.C*Calv) 

2  2 


[5-38] 


Due  to  the  complexity  of  the  coefficient  in  Equations  5-32  and  5-38,  these  equations  are 
solved  numerically  using  finite  differences,  liie  breath  is  divided  into  small  time  incremrats 
At  which  correspond  to  the  time  increments  used  in  the  numerical  calculation  of  the  regional 
flows  by  the  Ventilation  Model.  The  gas  concentrations  are  numerically  calculated  in  each 
r^ion  at  each  axial  position  at  the  end  of  each  time  increment. 

In  order  to  solve  the  differential  equation,  initial  conditions  and  boundary  conditions  are 
required.  The  initial  condition  is  simply  the  tracer  gas  concentration  in  the  lung  b^ore 
inhalation  starts.  At  the  distal  end  of  the  model  a  no-flux  boundary  condition  is  imposed  as 
no  tracer  gas  can  pass  through  the  back  wall  of  the  trumpet  channel  shown  in  Figure  5-2. 
During  inhalation,  the  concentration  of  the  tracer  gas  at  the  end  of  the  trachea  provides  the 
boundary  condition  at  the  entrance  to  the  channel,  but  this  concentration  is  not  known  and 
must  be  calculated  as  a  function  of  time.  In  order  to  do  this,  the  mouth  and  pharynx  are 
treated  as  an  extra  length  added  to  the  trachea  which  is  represented  as  a  rigid  tube  of 
constant  cross-sectional  area.  Equation  5-32  gives  the  mass  balance  in  this  tube  with  Jalv  = 
0.  This  differential  equation  is  solved  numerically  in  order  to  calculate  the  concoitration  of 
the  tracer  gas  at  the  end  of  the  trachea,  the  entrance  to  the  four  trumpet  charmels  that  form 
the  Convection-Diffusion  model. 

5.3.3  Results  and  Validation. 

Figure  5-17  shows  the  data  collected  during  a  single  breath  Nj  washout  test  (Nj  fractions  are 
dry  gas  fractions)  as  well  as  the  concentrations  calculated  at  the  mouth  during  exhalation  in 
the  simulation  performed  with  the  Pulmonary  Model.  In  a  single  breath  Nj  washout  test,  a 
subject,  starting  at  a  preset  lung  volume,  inhales  a  tidal  volume  of  gas  that  is  low  in 
usually  pure  Oj,  then  exhales  to  RV.  During  exhalation  the  N2  concoitration  at  the  mouth  is 
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measured  and  the  plot  of  N]  concentration  vs.  exhaled  volume  (or  time)  is  called  the  washout 
curve.  In  this  simulation,  the  measured  flows  were  used  as  an  input  in  the  Ventilatiim 
Model.  The  r^onal  flows  and  volumes  from  the  Ventilation  M(^el  were  input  into  the 
Convection-Diflfu5don  Model  and  the  N,  concentration  throughout  the  four  regions  was 
calculated.  The  subject  that  performed  this  single  breath  Nj  washout  test  had  an  RV  =  1600 
cm^  and  a  TCL  ==  7040  cm^. 

The  calculated  N2  washout  profile  compares  well  with  the  measured  N2  washout  curve.  The 
gas  originally  in  the  instrument  deadspace,  the  mouth,  pharynx  and  trachea  is  the  first  gas  to 
enter  the  lungs.  Once  this  deadspace  gas  has  been  inh^ed  fresh  tidal  gas  starts  to  enter  the 
lungs  and  a  front  is  quickly  established  between  the  inhaled  gas  which  has  a  low  N2 
concentration  and  the  resident  gas  which  has  a  high  N2  concentration.  Figure  5-18  shows  the 
calculated  concentration  profiles  at  1  sec  intervals  in  Region  1  as  a  function  of  the  summed 
volume  along  the  center  axis  of  the  channel.  The  regional  volume  and  the  inhaled  volume 
are  given  for  each  time  in  the  inset. 

As  inhalation  proceeds  the  N2  concentration  decreases  as  the  alveolar  gas  in  Region  1  is 
diluted  with  02,  but  the  center  of  the  front  that  separates  the  low  and  high  concentration  Nj 
gas  remains  near  50  cm^  (generation  15)  during  the  entire  inhalation  period.  The  axial 
position  of  the  front  is  determined  by  the  cross-sectional  area  of  the  trumpet  channel  and  the 
flow  rate.  The  front  is  at  the  position  where  the  N2  transport  by  convection  into  the  region  is 
equal  and  opposite  the  transport  by  diffusion  out  of  the  region,  i.e.  the  N2  is  convected  into 
the  region  as  fast  as  it  diffuses  mouthward.  In  the  airways  located  mouthward  of  generation 
16,  the  cross-sectional  area  of  the  channel  is  small  and  the  large  gas  velocities  in  these 
generations  makes  convection  the  dominant  transport  mechanism.  In  the  generations  located 
distal  to  the  front,  the  channel  has  a  very  large  cross-sectional  area  and  molecular  diffusion  is 
the  dominant  transport  mechanism. 

During  inhalation,  a  small  concentration  gradient  exists  in  the  first  5(X)  cm’  of  the  alveolar 
gas  in  Region  1,  but  the  gas  located  distal  to  this  position  is  completely  mixed.  Molecular 
diffusion  quickly  reduces  the  concentration  gradients  that  are  created  by  the  influx  of  into 
the  deep  alveolar  gas.  The  combination  of  the  large  cross-sectional  areas  and  the  short 
lengths  of  each  generation  enables  molecular  diffusion  to  be  so  effective.  The  distance  along 
the  center  axis  of  the  channel  between  a  volume  of  50  cm’  and  the  distal  end  of  the  channel 
is  only  1  cm. 

When  exhalation  starts  convection  no  longer  opposes  diffusion  and  both  mechanisms 
transport  the  N2  mouthward.  This  coupling  of  convection  and  diffusion  quickly  eliminates 
any  concentration  gradients  in  the  alveolar  gas  in  Region  1.  This  is  seen  in  Figure  5-18  at 
t=4  sec.  After  44  cm’  of  gas  have  been  exhaled  a  non-negligible  concentration  gradient 
exists  in  only  the  first  150  cm’  of  gas,  the  other  1170  cm’  of  gas  in  the  region  are  completely 
mixed.  At  t=5  sec  the  region  has  exhaled  only  127  chi’  of  gas,  but  all  of  the  remaining  gas 
is  completely  mixed,  there  are  no  concentration  gradients  within  the  region.  The  situation  in 
the  other  three  regions  is  similar,  except  the  N2  concentrations  are  lower  because  the 
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ventilation  to  volume  ratio  is  higher  in  the  more  dependent  regions  (see  Figure  5-10).  This 
means  that  more  gas  low  in  N2  enters  these  regions,  reducing  the  alveolar  Nj  ctxicentration 
below  those  seen  in  Region  1. 

While  the  agreement  between  the  measured  and  calculated  curve  in  Figure  5-18  is  good,  it 
turns  out  that  the  use  of  four  compartments  has  little  affect  on  the  washmit  because  the 
vertical  Nj  concentration  differences  are  relatively  small.  A  single  region  would  have 
produced  an  almost  similar  washout  curves,  but  the  regional  differences  are  important  in 
determining  the  alveolar  gas  concentration,  which  will  be  important  when  the  OKxiel  is 
expanded  to  include  transport  between  the  blood  and  alveolar  gas. 
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Figure  5-18.  Calculated  Nj  Fraction  in  Region  latlSec  Intervals  and  at 
the  End  of  Inhalation  {3.4  Sec).  X-Axis  is  Summed  Volumes 
Along  the  Center  Axis  of  the  Channel  that  Represents  Region  1 . 


77 


Regional  Volume 


SECTION  6.0 


79 


6.0  Cardiovascular  Model. 


6.1  General  Description. 

As  stated  previously,  the  cardiovascular  model  described  herein  parallels  the  developmoit 
presented  in  a  paper  by  White,  et  al.,‘^  as  well  as  papers  by  Jaron,  et  al.,'’-”'^*  and  earlier 
papers  by  Rideout,^  Snyder,^  and  Avula.^  Many  of  the  parameter  values  for  the  physical 
properties  of  the  vasculv  segments  were  taken  from  these  papers  as  well  as  Fung!*^  and 
Bergel.^  However,  some  of  the  parameters  listed  in  Jaron**  were  clearly  in  error  and  some 
of  those  employed  by  other  authors  were  not  useful  because  of  differences  in  the  layout  of 
their  models  and  ours.  Where  new  parameters  were  required,  they  were  derived  to  yield 
generally  acceptable  flow/pressure/volume  and  compliance  characteristics  of  the  various 
cardiovascular  subdivisions.  In  particular,  flow,  pressure,  and  compliance  data  from 
Burton,^’  and  Guyton,^*  were  employed. 

Figure  6-1  displays  the  general  schematic  for  the  model  compartments  and  their  connections. 

It  appears  feasible  and  desirable  to  extend  the  Phase  I  model  in  Phase  n  to  include  the  short 
term  physiologic  control  and  the  effects  of  pressure  breathing  on  cardiovascular  flows  and 
pressures.  This  will  enable  the  model  to  more  futhfiilly  simulate  the  true  response  of  the 
cardiovascular  system  in  response  to  short  term  changes  in  acceleration  and  breathing  rate, 
depth,  and  pressure.  An  enhanced  Phase  n  cardiovascular  model  could  be  validated  against 
experimental  data  and  then  employed  to  estimate  the  physiologic  effects  of  changes  in  the 
operating  conditions  of  the  breathing  system. 

6.2  System  Equations. 

The  model  describes  the  spatial  and  temporal  variation  in  the  mean  blood  pressure  along  the 
z-axis  of  the  body.  The  present  implementation  neglects  the  non-linear  and  convective  terms 
in  the  Navier-Stokes  Equation  as  well  as  the  complication  of  radial  flow  (flow  perpendicular 
to  the  vessel  wall;.  The  flow  is  assumed  to  be  laminar  except  in  the  ascending  an  descending 
aorta  where  fluid  flow  resistances  were  multiplied  by  33  to  account  for  increased  pressure 
losses  in  turbulent  flow.^*  The  following  set  of  coupled  differential  equations  were  solved 
for  each  arterial  and  venous  vascular  segment. 


^  ^  -  «»U  *  -  WU  [6-391 

^  j  -  *•«'»  [6^1 
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•«?(0^  -  Q(fW 


[Ml] 


drjt)  . 

dt 


2  *11  ‘KO  'i 


Where, 
t  =  Time 

r  =  Radius  of  vascular  segment 
1  =  Length  of  vascular  segment 

P  =  The  pressure  in  the  segment 
Q  =  The  s^mental  volume  flow 
C  =  The  capacitance  of  the  segment 
L  =  The  inertance  of  the  s^ment 
R  =  The  viscous  flow  resistance  in  the  segment 
R,  =  The  vessel  expansion  resistance 
Po,  =The  hydrostatic  pressure  diflfeience 
across  the  segment  because  of  gravity 


[sec]. 

[mj. 

[m]. 

[Pa]. 

[m’-sec’*]. 

[m’Pa-‘]. 

pcg-rn'^  or  [Pa-sec^*m’] 
[Pa*  sec- m’]. 
[Pa*sec*m'’]. 


[Pa]. 


And,  the  following  approximations  for  R,  L,  and  C  were  taken  from  Rideout,  et  al.* 


Ol-H-l 

[M2] 

8*11  -r* 

9-Ul‘ 

[6-43] 

4*K 

^  3‘r'V'l 

[6-44] 

2 -Eh 

^  C 

[6-45] 

Where, 


E  =  Young’s  modulus  for  vessel  wall  [Pa], 

h  =  Vessel  Wall  thickness  [ro]. 

Po  =  density  of  blood  p^-m  ’]. 

fiQ  =  viscosity  of  blood  [N*sec*m^. 
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Finally, 


\  =  pQ'G^'go'l'Cos{9) 


[6-46] 


Where, 

Gz  =  The  z-axis  “G-level”  in  units  of  earth’s  gravity  (unitless] 

go  =  The  earth’s  gravitational  acceleration  [m-sec'^ 

6  =  The  angle  between  the  s^moit  and  the  z-axis  [radians]. 

For  each  vascular  segment.  Equations  6-39,  6-40,  and  6-41  must  be  solved  simultaneously. 

Thme  are  20  artmial  and  20  venous  vascular  s^ments  (the  pulmonary  circuit  is  s^ment  1). 
There  are  10  peripheral  capillary  bed  s^ments  modelled  by  differential  equations  for 
pressure  and  flow.  The  periph^  capillary  resistance  and  oqncitaiioe  are  nooddled  as  ’’T” 
circuits  with  the  peripheral  capacitance  in  parallel  with  an  into  and  outlet  resistance  which 
represent  the  resistance  of  the  arteriole  and  venule  respectively.  The  pulmonary  circulation 
is  modelled  as  a  single  peripheral  bed  fed  by  a  single  arterial  segment  and  drai^  by  a 
single  venous  segment.  The  chambers  of  the  heart  are  modelled  as  variable  capacitances 
sq)arated  by  one-way  valves.  The  pulmonic  and  aortic  valves  are  also  modelled  as  cme-way 
valves.  The  pressure  reference  for  the  model  is  located  at  the  tricu^d  valve  which  is 
presumed  to  track  intrathoracic  pressure.  The  pressures  and  flows  in  the  various  s^ments 
are  coupled  by  their  spatial  connections.  The  z-axis  coordinates  for  locating  the  origin  and 
termination  of  each  segment  were  based  on  a  177  cm  tall  standing  man.  Including  the 
equations  for  pressure  and  flow  in  the  heart,  there  are  174  differential  equations  in  the 
model.  The  solution  was  computed  using  a  Runge-Kutta  numerical  integration  algorithm.*’ 

6.3  Simulation  Results. 

The  initial  conditions  for  segmoital  flow  were  set  based  on  a  steady-state  solution  for  the 
model  at  1  G.  (tr»(r<werse)  corresponding  to  a  supine  posture.  Initial  pressures  were  set  by 
computing  the  r  drc  r.zatic  pressure  at  each  segment  corresponding  to  the  initial  local  z-axis 
acceleration.  Po  c  is  adjusted  by  setting  the  appropriate  initial  orientation  angle  for  each 
segment  relative  to  the  z-axis  acceleration.  Postui^  and  acceleration  changes  during  a 
simulation  may  be  made  by  adjusting  the  orientation  angles  and/or  acceleration  during  the 
simulation,  lliis  can  be  most  easily  accomplished  by  adding  time  varying  profiles  for  G^  and 
the  B's  from  which  the  Pq^  values  for  the  flow  derivatives  are  calculated.  However,  the 
current  version  of  the  model  does  not  include  provisions  for  simulating  the  effects  of  the 
short  term  cardiovascular  reflexes  responsible  for  controlling  the  blood  pressure  and 
distribution  of  blood  flow.  Therefore,  no  attempt  was  made  to  match  transient  responses  of 
the  model  to  those  expected  in  the  response  of  the  true  cardiovascular  system.  Rather,  the 
model  was  coded  so  that  the  cardiac  output  was  maintained  at  a  constant  value 
(5.4  liters* min  ')  until  the  regional  blood  flow  values  returned  to  their  initial  values.  The 
resulting  steady-state  pressures  or  ’’equilibrium  pressures”  in  the  various  segments  represent 
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the  pressures  required  to  maintain  ‘‘normal”  blood  flow  through  the  vascular  and  peripheral 
segments.  The  results  of  the  two  simulations  described  bdow  illustrate  the  results  of  this 
process. 

To  illustrate  the  use  of  the  model  to  estimate  the  effects  of  acceleradcm  on  the  cardiovascular 
system,  three  simulations  were  conducted.  Figures  6-2  and  6-3  presmt  a  comparison  of  the 
equilibrium  arterial  and  venous  pressures  at  + 1  at  two  different  postures:  upright  and 
reclined.  The  segment  numbers  have  been  arranged  so  that  the  pulmonary  s^ment  (S^ment 
1)  is  at  the  far  left  and  the  remaining  segments  show  the  pressures  b^inning  at  the  most 
headward  segment  (Segment  20)  and  ending  at  the  most  footward  s^moit  (Segmoit  14). 

The  segment  orientations  for  the  reclined  posture  were  chosoi  to  iq>proximate  the  seated 
posture  in  a  high  performance  fighter  cockpit  with  a  seat  back  recliiiation  angle  of  20°  with 
the  feet  and  thighs  slightly  elevated  in  relation  to  the  buttocks.  The  second  set  of  simulaticm 
results  is  shown  in  Figures  6-4  and  6-S  which  illustrate  the  equilibrium  arterial  and  voious 
pressures  in  cardiovascular  system  for  the  reclined  posture  at  two  different  acceleration 
levels:  -i-1  Gz  and  -1-4  G^.  As  shown,  at  the  elevated  acceleration  the  more  headward 
arterial  segments  are  showing  negative  pressures. 
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Hgur«  6-3.  EquiHbritJin  Venous  Pressures 
Upright  +1  Gz  vs  Seated  -i-  1  Gz 


[6huiui]  untMUd  tnou«A 
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Figure  6-4.  EquiNbriuin  Arterial  Preeeures 
Seated  -i-l  Gz  versus  Seated  +4  Gz 


[Bhuhu]  iiintMid 


[«d]  Mnssajd  IVM^V 
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Figure  6-5.  Equilibrium  Venous  Preseures 
Seated  -i- 1  6z  versua  Seated  -f4  Gz 
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SECTION  7.0 
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7.0  Fhysiol<^  Data  and  Validation. 

7.1  Rapid  Decompression  Data. 

BRC  reviewed  the  rapid  decompression  data  furnished  by  the  Government.  This  series  of 
experiments  was  refmed  to  as  the  “EONS  Experiments”  because  they  were  carried  out  to 
verify  the  US  Navy’s  Expeditionary  Oxygen  Nitrogen  System  (EONS).  The  tests  were 
designed  to  study  ^e  effect  of  breathing  a  gas  mixture  with  less  than  the  98%  oxygen  at  high 
altitude.  The  motivation  for  testing  lower  oxygen  concentrations  than  found  in  Aviator’s 
Breathing  Oxygen  (ABO)  was  to  develop  evidence  that  gas  mixtures  produced  by  Molecular 
Sieve  Oxygen  Generation  Systems  (MSOGS)  were  acceptable  for  use  after  i^qnd 
decompression  to  S0,(X)0  feet  (15.24  km)  altitude.  The  studies  were  carried  (Hit  during  1986- 
1988  in  the  Altitude  Chambers  at  Brooks  AFB,  TX  by  the  staff  of  the  Crew  Technology 
Division  of  the  USAF  School  of  Aerospace  Medicine. 

BRC  obtained  the  raw  data  files  and  reviewed  them  for  suitability  as  data  for  validating  the 
ABS  m(xlel.  A  total  of  170  experimental  data  files  were  review^  and  assessed  for  quality. 
The  data  have  been  organized,  cataloged,  and  collated  by  subject,  by  experimental  condition, 
and  by  variable  for  easy  review  and  access  (See  Appendix  A).  The  Table  7-1  entitled, 
“EONS  Flights”  is  a  matrix  that  summarizes  all  of  the  test  data.  Each  row  of  the  table 
represents  an  individual  test  subject  and  each  column  represents  the  test  number.  Each 
element  of  the  matrix  contains  the  test  date,  the  analog  tape  number,  the  digitized  ASCII  file 
name,  and  notes  on  test  events.  Spreadsheet  files  have  b^n  created  in  Quattro  Pro  which 
summarize  the  test  data  according  to  test  subject  and  test  number.  Included  in  these  files  are 
time  history  plots  of  mask  pressure,  chamber  pressure,  partial  pressure  of  oxygen  and  carbon 
dioxide,  and  inspiratory  flow. 

The  test  subjects  wore  standard  USAF  fiight  gear  including  a  flight  suit,  and-G  suit,  helmet 
and  MBU-12P  oxygen  mask.  The  decompression  profile  consisted  of  a  five  minute  exposure 
to  20,(X)0  feet  altitude  while  breathing  either  dilut^  or  undiluted  supply  gas  from  a  standard 
USAF  CRU-73A  Diluter-Demand  Regulator.  At  the  end  of  the  five  minute  equilibration 
peritxi,  the  chamber  was  decompressed  to  50,0(X)  feet  in  approximately  1  sec.  The  chamber 
altitude  was  then  held  for  1  minute  before  a  descent  to  40,()()0  feet  was  initiated.  The  data 
presented  in  the  charts  shows  the  time  history  of  the  recorded  variables  for  a  two  minute 
period  beginning  one  minute  before  the  decompression  and  ending  one  minute  after  the 
decompression. 

The  respiratory  variables  recorded  in  the  data  files  were  inspiratory  flow,  mask  cavity 
pressure,  altitude  chamber  pressure,  and  the  partial  pressures  of  Oxygen  (PO2)  Carbon 
Dioxide  (PCOj)  at  the  end  of  each  breath  (end-tidal).  Thus,  the  data  are  not  (x>ntinuous,  but 
are  sampled  from  the  original  continuous  analog  tapes  at  the  end  of  each  breath.  The 
inspiratory  flow  data  were  supplied  as  the  peak  inspiratory  flow  from  the  breath  prior  to  the 

sampling  time  and  as  time  averaged  inspiratory  flows  from  previous  breaths  (  ).  The 
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time  averaging  was  achieved  with  a  three  breath  moving  average  filter.  The  inspintory  flow 
data  plotted  in  the  charts  are  the  time  averaged  data,  although  both  measures  were  included 
in  the  data  files  and  spreadsheets.  To  illustrate  the  character  of  the  data,  examples  of  data 
from  Subject  BES  are  shown  in  Figure  7-1,  below.  In  the  graph,  the  abscissa  is  time  in 
seconds,  the  left  ordinate  is  absolute  pressure  and  the  right  ordinate  is  flow.  The  chart  title 
encodes  the  experimental  conditions.  In  Figure  7-1,  the  data  are  displayed  for  subject  BES, 
breathing  undiluted  100%  oxygen  before  and  after  the  5  psig  decompression  from  20,000  to 
50,000  feet  (6096  to  15240  m).  The  data  show  the  chamber  pressure  and  mask  pressure 
falling  immediately  after  the  decompression  begins.  The  chamber  pressure  fells  to  88  mm Hg 
while  the  mask  pressure  falls  only  to  slightly  over  1 10  mmHg  because  of  the  onset  of  30 
mmHg  pressure  breathing  (PBA).  The  end-tidal  PO2  falls  from  approximately  260  mmHg 
before  the  decompression  to  50  mmHg  post  decompression  while  the  corresponding  PCOj’s 
fell  from  approximately  35  mmHg  to  30  mmHg.  The  inspiratory  flow  is  seen  to  have  risen 
briefly  from  about  7  liter- min  *  after  the  decompression  and  then  return  to  the  15  to  20 
liter- min  *  toward  the  end  of  the  exposure. 
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Figure  7-1.  Subject:  BES/100%  02  Non-Dilution 
20/50  kft  Rapid  Decompression 
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RD  =  Rapid  Decompression,  ML  =  Mask  Leak,  MS  =  Mask  Spec,  LOC  =  Loss  of  Consciousness 
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RD  =  Rapid  Decompression,  ML  =  Mask  Leak,  MS  =  Mask  Spec,  LOC  =  Loss  of  Consciousness 
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RD  =  Rapid  Decompression,  ML  =  Mask  Leak,  MS  =  Mask  Spec,  LOG  =  Loss  of  Consciousness 


7.2  Use  of  EONS  Data  for  Model  Validation. 


The  various  modules  which  will  eventually  comprise  the  integrated  ABS  model  have  been 
individually  run  under  representative  conditions  to  demonstrate  their  agreement  with 
physiologic  data  from  die  literature.  However,  the  ABS  model  cannot  be  fully  validated 
until  the  component  modules  are  integrated. 

The  EONS  data  presented  in  Appendix  A  will  be  essential  in  validation  of  the  ability  of  the 
ABS  to  predict  the  transient  flow,  pressure,  and  respiratory  gas  composition.  Bomar^  has 
previously  demonstrated  the  use  of  the  EONS  data  to  vali^te  a  simple  dynamic  model  of  the 
changes  in  partial  pressure  of  respiratory  gases.  Figure  7-2  show  the  end-ddal  oxygen  and 
carbon  dioxide  partial  pressures  (PejOj  and  PetCOj)  predicted  by  his  model  (dashed 
continuous  lines)  and  those  which  occurred  during  experimental  rapid  decompression  of 
human  subjects  breathing  100%  and  93%  oxygen  diluted  with  air  (symbols).  There  was 
reasonable  agreement  between  predictions  and  the  data  which  improved  with  post¬ 
decompression  elapsed  time.  This  was  true  despite  the  fact  that  Bomar  did  not  model 
lung/chest  mechanics  and  he  assumed  constant  rates  of  oxygen  consumption  and  carbon 
dioxide  production  (i.e.  O2  uptake  and  CO2  excretion  were  assumed  to  be  independent  of  the 
partial  pressure  gradients  across  the  alveolar  membrane).  These  initial  results  were 
encouraging  and  pointed  to  specific  modifications  required  to  bring  the  ABS  model  in  closer 
agreement  between  the  predictions  and  the  experimental  data.  The  Phase  II  model  will 
include  partial  pressure  dependent  alveolar  gas  transport  as  well  as  a  dynamic  cardiovascular 
model  which  will  mimic  the  transient  changes  in  pulmonary  blood  flow  as  well  as  the 
distribution  of  b.c.-!  flow  along  the  vertical  axis  of  the  lung.  This  will  enable  ventilation 
perfusion  ratio  estimates  to  be  made  for  comparison  with  data  from  the 
literature.*’^*  ■*^’*^’^’^*’** 

The  EONS  data  will  also  allow  a  limited  validation  of  the  model  for  predicting  the 
physiologic  effects  of  PBA  at  30  mmHg.  Other  than  the  EONS  series,  there  are  data  from 
Emsting”  on  the  major  effects  of  pressure  breathing  on  the  respiratory  and  cardiovascular 
systems.  On  the  other  hand,  there  are  only  limited  data  on  the  effects  of  pressure  breathing 
on  ventilation-perfusion  relationships.**  *’  Undoubtedly,  the  attempt  to  validate  the  ABS 
for  higher  breathing  pressures  and  its  effect  on  ventilation-perfusion  will  point  to  the  need  for 
specific  data  so  that  experimental  protocols  can  be  suggested. 
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Figure  7-2.  Rapid  Decompression  20/50  kft  100%D  02 
End-Tidal  02  and  C02  Partial  Pressures 


nme  (sec  [RD@(=#)) 
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8.0  Conchisioiis  and  Recommendations. 

8.1  Conclusions 

BRC  completed  an  SBIR  Phase  I  project  to  study  the  feasibility  of  developing  a  model  of  the 
Aviator’s  Breathing  System  (ABS).  The  motivation  for  the  project  was  the  desire  to  develop 
a  model  which  could  simulate  the  cardiovascular  and  respiratory  responses  to  altitude  and 
acceleration  stress  encountered  in  high  performance  military  aircraft.  Once  validated,  the 
model  could  be  employed  to  integrate  existing  physiologic  knowledge  and  to  estimate  the 
physiologic  effects  of  changes  in  the  environment  or  the  design  of  aircrew  protective 
equipment. 

The  Phase  I  effort  was  concentrated  in  developing  computer  software  modules  which  were 
capable  of  simulating  the  main  effects  of  environmental  variables  on  the  breathing  system 
hardware  and  respiratory  and  cardiovascular  systems.  Software  modules  were  developed  and 
tested  for  simulation  of:  (1)  the  flows  and  pressures  within  the  breathing  gas  delivery 
system;  (2)  the  flows,  pressures,  and  respiratory  gas  distribution  within  the  lung;  and  (3)  the 
steady-state  flows  and  pressures  within  the  cardiovascular  system.  Subprograms  were  also 
developed  to  compute  altitude-barometric  pressure  relationships  as  well  as  passenger  cabin 
pressures  in  military  aircraft. 

In  addition  to  the  software  development,  BRC  reviewed  and  organized  the  Government 
furnished  data  from  a  series  of  manned  rapid  decompressions  known  as  the  “EONS 
Experiments."  The  data  from  approximately  170  experimental  decompressions  were 
screened  for  their  suitability  for  use  in  parameter  selection  and  validation  of  Lhe  .  ..spiratory 
modelling  software.  The  data  appears  to  be  highly  coherent  and  fully  usable  for  model 
validation.  The  summarized  results  of  the  EONS  Experiments  represent  a  new  physiologic 
database  which  is  directly  applicable  to  USAF  operational  scenarios.  The  respiratory  gas 
composition  time  histories  represent  original  scientific  data  suitable  for  publication  and  use 
by  the  aerospace  medicine  community  independently  of  their  value  in  the  ABS  modelling 
effort. 

The  work  accomplished  in  Phase  I  has  laid  the  groundwork  for  developing  an  integrated  ABS 
Model  in  a  Phase  II  SBIR  program.  BRC  has  written  prototype  software  modules  which 
form  the  major  building  blocks  for  the  Phase  II  ABS  Model.  We  conclude  that  the 
development  of  an  integrated  ABS  Model  is  feasible  and  desirable.  We  also  conclude  that 
there  is  sufficient  existing  data  for  selection  of  parameters  and  independent  validation  of  the 
major  components  of  the  Phase  II  ABS  model.  However,  the  complete  validation  of  the 
combined  effects  of  acceleration,  pressure  breathing  and  altitude  on  the  ABS  must  await  new 
experimental  data.  Fortunately,  the  attempt  to  simulate  the  physiologic  response  to  combined 
stressors  will  explicitly  point  to  missing  data  and  guide  the  development  of  experimental 
protocols  to  elucidate  their  significance.  The  following  section  outlines  the  major  steps 
required  to  create  the  Phase  II  model  by  building  on  the  work  accomplished  in  Phase  1. 


8.2  Recommendatioiis. 


At  present,  none  of  the  Phase  I  ABS  modules  is  fully  validated.  Therefore,  prior  to  their 
integration,  additional  development  of  the  Phase  I  software  modules  will  be  required  to 
ensure  their  ability  to  independently  simulate  the  physiologic  response  to  changes  in  their 
input  variables.  In  addition,  the  software  code  and  numerical  methods  for  solution  of  the 
system  equations  will  be  optimized  and,  where  possible,  standardized  across  modules.  This 
will  result  in  more  efficient  code  and  ease  the  process  of  creating  the  software  interfaces 
necessary  to  integrate  the  Phase  I  modules  into  a  fully  integrated  ABS  Model.  The  following 
sections  outline  the  major  work  on  each  software  module  required  before  their  integration. 

8.2.1  Breathing  Gas  Delivery  System. 

The  breathing  gas  delivery  system  is  comprised  of  a  source  of  breathing  gas,  such  as  an 
oxygen  supply  or  a  Molecular  Sieve  Oxygen  Concentrator  (MSOC),  a  demand  regulator 
which  regulates  the  flow  and  pressure  of  the  breathing  gas  in  response  to  the  changes  in  the 
ambient  pressure  and  acceleration,  and  finally  an  oronasal  mask  which  ensures  unidirectional 
flow  of  breathing  gas  and  forms  a  seal  to  the  face  of  the  aviator. 

The  preferred  way  to  model  the  breathing  gas  source  would  be  to  model  the  composition  at 
the  demand  regulator  inlet  as  a  function  of  the  mass  flow  drawn  by  the  regulator  and  the 
ambient  pressures.  This  approach  will  work  for  both  MSOC  based  systems  or  systems  based 
on  storage  of  gaseous  or  liquid  oxygen.  Backup  oxygen  systems  such  as  seat  mounted 
Emergency  Oxygen  Systems  are  relatively  simple  to  model.  Regulated  Back-Up  Oxygen 
Systems,  such  as  those  employed  with  MSOCS  can  be  modelled  similarly  to  the  demand 
regulator,  but  without  dilution  of  the  feed  gas. 

The  demand  breathing  gas  regulator  in  the  present  model  assumes  the  flow  demanded  is 
supplied  at  a  pressure  typical  of  the  functional  performance  of  an  ASCC  demand  regulator. 
For  the  Phase  II  regulator  model,  it  will  be  necessary  to  completely  parameterize  the  mass 
flow-pressure  characteristics  of  the  regulator  and  include  G-sensitivity  as  required. 
Parameterization  of  the  regulator’s  performance  is  preferable  to  any  attempt  to  model  the 
internal  hardware  function.  Breathing  gas  regulators  are  highly  nonlinear  devices  and  the 
eflfort  required  to  write  and  solve  the  equations  describing  the  flow-pressure  characteristics  of 
even  a  simple  regulator  would  be  substantial.  The  regulator  model  will  also  include  a 
parametric  representation  the  regulator’s  supply  gas  dilution  characteristics  as  function  of 
ambient  and  outlet  pressures.  The  regulator  outlet  relief  valve  will  be  modelled  on  the 
characteristics  of  a  poppit  relief  valve  rather  than  the  ideal  relief  included  in  the  present 
model.  The  mask  tube  and  connector  flow-pressure  characteristics  will  employ  the 
isothermal  tube  flow  equations  and  will  allow  for  turbulent  pressure  losses. 

The  oronasal  mask  model  is  essentially  complete  as  implemented.  The  model  is  based  on 
the  functional  performance  of  an  ASCC  compliant  breathing  mask.  The  present  model  can  be 
modified  to  change  the  mask  valve  area  versus  pressure  function  to  match  an  actual  mask  if 
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desired.  Empirical  data  can  be  employed  to  match  the  valve  characteristics  in  the  model  to 
actual  mask  flow-pressure  data. 

8.2.2  Pulmonary  Model. 

The  ventilation  pattern  in  the  Pulmonary  Model  is  presently  driven  by  the  flows  at  the  mouth 
(see  Figure  5-1).  In  order  to  integrate  this  model  with  the  mask/regulator  model,  it  will  be 
necessary  to  make  pleural  and  ambient  pressures  drive  the  gas  flows  in  the  model,  i.e.  the 
pressure  difference  across  the  chest  wall.  This  will  require  adding  the  elastic  characteristics 
of  the  chest  wall  to  the  Pulmonary  Model  and  changing  the  logic  of  the  numerical  solution 
scheme.  This  addition  will  make  the  model  more  realistic  as  pleural  pressure  changes 
normally  drive  gas  flow,  and  in  the  case  of  a  rapid  decompression,  changes  in  the  ambient 
pressure,  coupled  with  the  pleural  pressure,  creates  gas  flow  in  the  lungs. 

Presently  the  model  assumes  that  the  ambient  pressure  remains  1  atm  and  that  the  density  of 
the  gas  in  the  lungs  remains  constant.  Obviously,  this  is  not  the  case  when  pressure 
breathing  and  rapid  decompressions  are  modeled,  so  the  Pulmonary  Model  will  require 
modification  to  include  the  density  effects  of  lung  gases  which  are  not  presently  included. 
This  can  be  done  with  relatively  little  difiSculty. 

The  Pulmonary  Model  inherently  contains  the  effect  of  acceleration  on  ventilation  through 
the  G'term  in  Equation  5-18,  but  the  resulting  regional  volumes  and  flow  patterns  for  G- 
values  greater  than  + 1  G*  have  not  been  studied  in  detail,  and  the  performance  of  the  model 
in  predicting  ventilation  patterns  in  high-G  environments  will  require  verification. 

In  order  to  integrate  the  Pulmonary  Model  with  the  Cardiovascular  Model,  the  transport  of 
carbon  dioxide  and  oxygen  between  the  alveolar  gas  and  the  blood  must  be  included.  The 
model  presently  tracks  the  concentration  of  a  single  gas  in  the  airways  of  the  four  regions, 
therefore  it  cannot  realistically  include  the  effects  of  more  than  one  tracer  gas  species.  To 
accommodate  alveolar  exchange,  the  model  must  be  modified  tr  '<ce  multiple  gas  species 
and  to  include  the  transport  of  oxygen,  nitrogen,  and  carbon  du  ;  throughout  the  airways 
and  alveoli  of  each  region  in  the  model. 

8.2.3  Cardiovascular  Model. 

The  major  change  required  in  the  cardiovascular  model  will  be  to  incorporate  dynamic 
response  characteristics  and  add  models  of  the  short  term  physiologic  control  systems.  The 
steady  state  effects  of  acceleration  are  already  implemented  in  the  cardiovascular,  but  the 
modifications  will  be  required  to  include  transient  changes  in  flow  pressure  within  the 
segments.  Pressure  breathing  effects  will  be  implemented  by  modelling  the  rise  in 
intrathoracic  pressure  and  its  effect  on  venous  return  and  blood  pressure.  Methods  of 
modelling  pressure  garments  such  as  partial  pressure  garments  and  anti-G  suits  must  also  be 
added.  It  should  be  feasible  to  model  the  application  of  external  pressure  by  manipulation  of 
the  tissue  pressure  (Pe«)  which  is  already  included  in  the  model.  The  current  version  of  the 
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model  also  runs  very  slowly  with  the  Runge-Kutta  numerical  integration  algorithm  presently 
employed.  Different  methods  of  formulation  and  solution  of  the  differential  equations  should 
be  investigated  to  optimize  the  speed  and  stability  of  the  model. 

The  pulmonary  cardiovascular  section  will  be  split  into  four  subsections  each  supplying  the 
corresponding  region  of  the  pulmonary  model.  Local  regulation  will  be  simulate  as  will 
flow  pressure  alterations  in  response  to  acceleration  and  pressure  breathing. 

8.2.4  Integration  of  Subsystem  Modules. 

Phase  1  has  created  basic  modules  and  a  framework  for  an  integrated  ABS  model.  As 
modem  life  support  systems  (LSSs)  are  substantially  integrated,  the  Phase  II  model  will 
attempt  to  include  the  interaction  of  the  relevant  LSS  components  with  the  environment  and 
the  crew  member  and  other  parts  of  the  LSS.  Most  of  these  variables  will  simulated 
parametrically  as  inputs  to  the  system  to  model  PBG  as  well  as  anti-G  suit  inflation  during 
pressure  breathing  for  altitude  (PBA).  Modelling  of  interfaces  will  likely  be  one  of  the  most 
interesting  and  most  difficult  parts  of  combining  the  Phase  I  modules  into  an  integrated  ABS 
model  in  Phase  11.  To  our  knowledge,  no  one  has  modelled  the  interfeces  between  the  life 
support  equipment  and  the  aviator.  Just  as  importantly,  with  the  aid  of  the  ABS  Model  it 
should  be  possible  to  discover  the  nature  of  improvements  which  would  be  necessary  to 
correct  poor  interaction  or  poor  integration  between  the  crew  member  and  the  LSS.  The 
integrate  ABS  Model  will  be  employed  to  perform  parametric  studies  as  an  aid  to 
illustrating  and  understanding  the  nature  of  the  interactions  between  the  environment,  the 
LSS,  and  the  aviator.  A  knowledge  of  the  (dynamic)  nature  of  the  interaction  between  the 
physical  and  physiologic  subsystems  will  also  help  guide  the  Government’s  research  program 
toward  elucidating  the  physiology  necessary  to  describe  the  effects  of  LSS  design  on  aircrew 
protection. 

Once  the  individual  modules  are  created  and  validated,  they  will  be  integrated  into  a  single 
software  program  or  an  interactive  suite  of  programs.  The  validation  of  the  integrated 
software  will  require  more  than  a  recomparison  of  simulations  to  the  same  data  employed  in 
the  validation  of  the  individual  models.  Data  will  be  required  which  describes  the  effects  of 
a  single  environmental  variable  on  both  the  pulmonary  and  cardiovascular  systems.  In 
addition,  data  on  the  response  of  the  pulmonary  and  cardiovascular  systems  to  combined 
environmental  inputs  will  be  required.  This  phase  of  die  model  validation  is  likely  to  be  the 
most  difiicult  because  of  the  paucity  of  data  on  such  effects.  Undoubtedly,  the  attempt  to 
validate  the  integrated  model  will  point  to  the  need  for  specific  new  data. 

The  design  of  user  and  software  modules  will  be  a  major  task  of  the  Phase  II  development 
program.  The  goal  will  be  to  design  a  single  software  module  with  a  “user  friendly” 
interface  which  allows  easy  selection  and  adjustment  of  input  variables  and  display  of 
simulation  results.  It  is  likely  that  the  display  of  model  time  histories  will  require  use  of 
separate  graphics  programs  or  spreadsheet  models. 
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Finally,  it  is  recommended  that  the  Phase  II  development  proceed  methodically  and  be  spread 
over  a  sufficient  time  to  allow  for  development  of  new  data  which  may  be  required.  The 
Phase  II  models  should  be  developed  as  independent  modules  for  use  independently  of  the 
integrated  ABS  model.  This  will  permit  “stand  alone”  use  of  the  modules  or  for  ffieir 
further  development  for  projects  outside  the  context  of  the  ABS. 
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